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SUMMARY 


The problem of dense output (luring the numerical solution of an ordinary 
differential equation, using a Hunge-Kutta algorithm, oan severely reduce the 
effioienoy of the integration method* A set of soaled Runge-Kutta algorithms 
for the third- through fifth-orders are developed to determine the solution at 
any point within tho integration step at a relatively small lnoreaae in com- 
puting time* Baoh soaled algorithm is designed to be used with an existing 
Rungo-Kutta formula, using the derivative evaluations of the defining algorithm 
along with an additional derivative evaluation (or. two)* Third-order, soaled 
algorithms are embedded within the existing formulas at no additional deriva- 
tive expense* Such algorithms oan be easily adapted to generate interpolating 
polynomials (or dependent variable stops) efficiently. 


INTRODUCTION 


The Runge-Kutta algorithm, designed to treat the initial value problem 


^ * f(t,y) , y(t 0 ) * y 0 

at 


(D 


where y is an n-veotor, is an explicit single-step algorithm. A linear combi- 
nation of derivative evaluations to approximate a trunoated Taylor series is 
used. These algorithms are frequently used in solving the ordinary differen- 
tial equation (ODE) because of their high aoouraoy and simplicity. For Runge- 
Kutta methods to perform efficiently, however, the algorithms must be operated 
at an optimal step size that is consistent with the truncation error oriteria. 

A significant reduction in effioienoy may be encountered if frequent data 
output is requested. This difficulty arises because the Runge-Kutta algorithm 
determines the solution only at the end points of the integration stsps; thus, 
the step size must be reduoed to ooinoide with the requested output points. 
Shampine, Watts, and Davenport (ref. 1) have shown, however, that reduction of 
the step size by this procedure seriously impairs the effioienoy of the method; 
i.e., the effioienoy of the Runge-Kutta algorithm depends on the ability to use 
large step sizes whenever permissible during the solution of the ODE. (Unlike 
the multistep methods, a backlog of solution points for an interpolating polyno- 
mial is generally unsatisfactory* This is due to the large step sizes used by 
the Runge-Kutta methods, which provide too wide a spacing for an aoourate 
interpolating polynomial.) 

An alternate prooedure for data determination may be developed. This pro- 
cedure allows the Runge-Kutta algorithm to be operated by using a near-optimal 
step size While still providing the solution at an intermediate point for only 
a slight increase in expense. Such a prooedure, first developed for Runge- 
Kutta-Nystroem methods (ref. 2), involves the development of a oompanion 
algorithm of equivalent order, using the funotion evaluations of the basic 
algorithm plus an additional evaluation (or two). The new scaled algorithms 
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givo the solution veotor at an Intermediate point t * t 0 ♦ h # where h # « oh 
set of ooeffioinnf OS ! tiVO p flCali S g pnPamoter usually loss than unity. pSr each 

«-jrjr.a®a 

of the KSrtSS T/oZZT ° f ooofflol ™ ta ’ » <*• oflloiMoy foaturoa 


»n.,x,4 The d ? rivafc i on of the soaled nlgorlthma and a list of ooeffloients for 
varying values of o Qro presented for tho embedded pairs of Runae-Kutta 

» ^hl^rrCrlf^SraSd^another 

have £, ^ z SS^^tTJS 

feL 15 ® ubr ^ uti ^ e r KEQ^ whioh evaluates truncation error coefficients through 
the 12 th-order for Runge-Kutta methods. These scaled ^Runge-Kutta oSeffio «n?^ 

Droo e eL Q iL C :S^ d dl / ectly f ™ the Senerating progJaL into a large cofcpSje 
processing system for reproduction. Similarly, the scaled coefficients listed 
in rational form hav » been computed with the aid of a rational arithmetic ni Lie 
age and verified by subroutine RKEQ before being copied into thl ZSS ItLlt 

SSJff-- # dls — 10 " *• implementation XI 


THE SCALED RUNGE-KUTTA ALGORITHM 


The basio Runge-Kutta algorithm of order 
U) assumes the form 


p for the solution of equation 


y (t 0 + h) a y 1 + 0 (hP +1 ) S y 0 + h ! C K f K + 0 (hP +1 ) 

KaO 


(2a) 


where 


fo a f(t o ,y 0 ) 


(2b) 


and 

« 


(' 


K-1 


f < 8 r t 0 + <*h, y 0 + h r p K x f x 

V XsO ’ 


(2c) 


Private communication. 

*xr Uon Error C08m ° 1 '” ts 
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- or pr-srs }•. -/*»?-- 

-ssifs ssjsjsts, to B ^ 

Z%o^^ 

higher order.) The equetlooe of oondltlon ere develop™ wHhlhe LSmonel 

» 


assumptions that 


K-1 

x »o e< * x ' “* «> 


for ic * 1,2,... # p, whioh define the B Kf0 ooeffioients. 

defining ° f the ® oaled algorithm requires knowledge of the 

r®; 1 " 111 ® algorithm because the new methods use the basio ooeffioients and 

a?L^^ at }r e , eValUatl<>n8 as the «*• of ^eir own sya?em? In thf soared 
U?e?-f? rl iS aa ? uae h tha J an aooeptable step of length h has been token 
min^d'at 1 '/!’***! ft ha !® bee " 1 evaluated and the solution has been deter- 
SS; “Jhe^e ° % ^ n ° W tha S ° lufcion at fc 8 to * h 8 should be detor- 


h* s a h 


(4) 


with o as some soaling faotor where, generally, a e (0 i) The MU - rta1a j 
wlth^he^solutlon^and ** *“*• 


* • P # » # 

yi 8 y 0 + h r qc f K (5a) 

•CsO 


where 


f o 8 f<Wo> (5b) 


and 
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• / • * * K-1 • « 

8 t l fc o ♦ °Kh , y 0 h £ f* 


0 Mi 

for K b l,3,...,p , If tha a and p* coefficients can be chosen auoh 
that 


f K 8 f K 


f?'Vj’ 2 *hen tho expense of applying the soaled algorithm will be 

minimal because fi, eg, ... , f p have already been oomputed. Comparing equa- 
tions (2b) and (2c) to equations P (5b) and (5c) show that 

« i 


f K 8 f K 


i i 

<*k h 


8|C, X h 8 8 k, X h 


• § 

for k = 1,2,...,p; X * 0,1 ,..., k- 1. Thus, only f D+1 , ... , f 6 » need to be 
oomputed. (While a third-order algorithm exists, requiring no additional de- 
rivative evaluations, a more aoourate third-order formula may be derived by 
using one additional evaluation. For the fourth-order algorithms, one addi- 

13 needed <P* 8 P ♦ D- For the fifth-order algorithms, two 
extra evaluations are necessary (p» a p + 2) for a fifth-order solution, al- 
though a fourth-order solution may be obtained for only one additional 
evaluation.) But because h« a a h, 


«K a OtK / o 


8k, X 8 8 k ,x / o 


f . y jo* ** 9 



for k *» 1,2,...,p} X ■ Thus, many of the coefflolentn for each 

scaled algorithm are already determined by being related to their corresponding 
ooeffloionta in the defining algorithm by a factor of 1/a. 

An in the development of Runge-Kutta algorithms, the equations of condi- 
tion for tho sealed methods may bo solved with one or more values of the 
<** or 0* coefficients an free parameters, Tho main criterion for selecting 
those parameters in tho minimisation of the norm of tho error oooffiolonts for 
tho first neglected term in tho Taylor sorios oxpannlon. More specifically, 
tho oocuracy of an algorithm of ordor p may bo Imprevod by reduoing tho error 
ooeffie lento for order p+1 so that 

«. f*p+1 «2 

°p+1 2 ^ Z T p+1,J 

l J= ’ 

approaohes a minimum, where Xp+T is the number of error ooeffioients of 
order p+1. As in the Runge-Kutta algorithm development, however, aeoondary 
conditions suoh as moderate bounds on. the ooeffioients, show that near-optimal 

• 

values of G p+ i may produoe algorithms that are more effioient than the 
algorithm for which G p+ i minimum is aohieved. Thus, in selecting the free 
parameters, algorithms having a value of Gp+i olose to G p+ i minimum are 
studied, but it is not demanded that G p+ i minimum be reached. 

The stability properties of the Runge-Kutta algorithms are also of oon- 
oern in the development of the ooeffioients. The stability analysis, applied 
to the single, linear differential equation 


1/2 p+1 

0 


(9) 


y s X y 


determines the largest value of the step size h for whioh errors in the solu- 
tion are not amplified by the algorithm itself, where X is any eigenvalue of 

the Jaoobian matrix of equation (1) with X oomplex. For an algorithm 

fiy 

of order p and a given X , the step size h is restricted by the polynomial 


P(Xh) 


P 


a E 
k=0 



k! 


+ 


p # +1 

E 


kep+1 


Pk (Xh) 


k 


( 10 ) 
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where depond a upon the 




0 *. 


f 2 ; |P<Xh)| i i: nd 0ftr?S 1 ^^^ f ;rrthe POr 

imagine™ aSt^ 9 * fl J ab J ut y °^vea oroaa the negative real nxla R^ and thn 
frJ a 2n?«rL ftXt0 1 V, In dmfQ l°P J nR the aoalod algorithms, oaro should bo taken 
advanced Bfc ? blU iH bQandfl i although, if the solution la not 
fehfl ^LJ ? fch ?.C nfc , fc ° * h » tba amplification during the atop taken by 

XJfcattons! 80 ftff ° Qfc ° nXy fch ° lmmodlat « °«tput «"d not the further 


THIRD-ORDER SCALED ALGORITHMS FOR RK<2)3 FORMULAS 


thoae C for f the°BKT(2)3 ^° T ‘ 3) ,p “ found ln tabl ° 1 “ hll » 

listed In (.Ihif TT *u 8orlfcbra havln « lower truncation error ooeffiolonta, are 

IZfiZL? "f22i; “: K **» l»iwd formulas my be uaed "ith either order. The 
fhi U iliHi.f 0rB ? Wi ^ h !' he C °9 0ffioi ®nts is of the lower (second) order and 
diffW™l°h f ormed " itb thQ ^ ooeffiolonta is of higher (third) order! (The 
h b J p" tW ° “^tlo** 8 Sives an estimate of the looal error that 
JSJ ‘ e h L^ 8 »^J?: f t0p ; alz ! » f ‘he lining algoritS. Soau^ 1 

mina/4 „k4 u 4 l fl a l*' 0ad y available, only one scaled solution is deter- 

eialust^n , 1 ? ° f .J he 5^ order.) By inserting one additional derivative 
of th^thirt o^rViJ^r’ SOalS ? ® lgorithm my be generated (although errors 
defining SLTSlf! l J beeause of the ^ion obtained by the 

tion ftnti m ^ L thm * ) An alternate solution, presented in the sec- 

algor ithm^at^nn^ona^ T !? i r d “° rder Solutions" generates a third-order, scaled 

nS» * at additional derivative expense. Beoause of the lack of free 

henee et fch?s ~ m f ybe zeroised over the truncation error coefficients: 

method Pon^M ed al6 ° ritb ® not be a particularly suitable third-order 
tive^valuatio^is^esented? re < Uiri " 6 °" e additlonal 

thlrd^d^^aScmhVSrn^l^ 1 " 8 the eqUatl °” S ° f «» 


» * * *2 

P 2,1 s 02,1 Ol *02/2 (11a) 



» • « « » 2 

03 ,iai ♦ 03 , 2<*2 *03/2 


(11b) 


and 


* * * • # § • §2 

p 4,1 * 04,101 ♦ 04,202 ♦ 34,303 * 04 f 2 


(11c) 
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# 

along with C 1 b 0 , (During the development of the defining algorithm, equa- 
felons (11a) and (11h) (without *) were imposed and, thus, are valid in the 
aoaled system, otherwise, sueh assumptions oould not bo arbitrarily imposed 
because the a" and eooffieients involved are already determined.) 


The paramofeorn for eonsferuoting the sealed algorithms are <14, ^4 i, 

i • ' 

&j| 3, and 84,35 that is, feho remaining oonffieiento aro either sealed 
values of feho defining oeoffioionta or, an in feho ease of feho tr eooffiolonts, 
are determined from these parameters. The C* eooffJelenfes may be writfcon an 


• «» «• • • i « • 

C2 S ( 1/14 - («3 + «n )/3 + «3 04/2) / (a 2 (a 2 - ot 3 )(a 2 - 014)) (1?) 


« 

with O3 


« 

and C4 determined by permuting indioes, and 


« 

Cq a 


§ 

The 84 coefficients are 


« • i 

-(03+03+04) 


84,1 a -< c 2 ^2,1 + C3 03 ,i)/C 4 


« 

® 4,2 8 


? 4,2 


§ §2 « »2 • • 

34 , 1*1 - 03(04 /2 - 64,1*1) 

<x*(o* - ot* ) 

2' 2 3' 


( 13 ) 


( 14 a) 


( 14 b) 


and 


« 

B 4,3 


s 


* » §2 a »2 * * 

P4,2 - 34,1*1 “ «2(«4 /2 “ 04,1*1) 


a*(a # - a # ) 
3 3 2 


( 14 c) 


7 


where 


• • * •?. « « §j? „ *? . 

P V * 0/12 ~ c 2 P 2,1 ft 1 - C3C»3,1 ft 1 ♦ >> / Ci, 


< 1G) 


W 

hf'Jh 0 !'° °?T fc 0 , V r ? e{ l uoUon (3). Ail other pernriietere ore determined 
by the gnnorafclng net of aaoffioionfca uolng oquntlonn (fln) nnd (fib). The aolu- 

tlon of the equations of condition in thon achieved with «! romnlnlng ho n 

pa JJ ,w i 0 J? th0 nmUd algorithm for the RK(2)3 In actually^? fourth 

ordor, tho defining algorithm ia of third-ordor aoouraoy (or nooond, if the 

aoourato^nn ?? lufclon J a used), and hence, the scaled oolutton will only bo an 
aeourato ao tho oompufcod solution at t a t 0l but at least additional orrora 

duo to the scaled algorithm will be of higher order. Because oj & o, the 

value of a 4 la selected to reduoe Gg (eq. ( 9 )), involving fifth-order error 

i 

ooeffijients Tg ( j, and to ensure reasonable magnitudes on all <**, 0 * 

and C ooeft loients (as far as possible). These error coefficients are rela- 

tively Insensitive to changes in <*4 f or both the rkf( 2)3 a nd the RKT(2)3 

formulas. Considering the RKT< 2)3 formula, a near minimal value of 0 * is 

achieved by selecting <*4 * 0.60 for 0 e ( 0 . 00 , 0 . 45 ], or . 0.50 for 

0 e ( 0 . 1 * 5 , 0.70], 04 = 0.30 for o t (0.70, 0.85], or <*J = 0.10 for 

• t 

a e (0.85,1.00). While O5 is nearly minimal for <*4 = 0.60 for 

V °»65] or <*4 S 0.50 for o e (0.65, 1 . 00 ) for the RKF( 2)3 scaled 

“ 3 ° aled w f olu J io « with aoouraoy equivalent to that of the defining 
be aohi ®ved at no additional derivative expense. (See seotion 
one^adfH t H ^® badd ® d Third-Order Solutions".) Therefore, the scaled RKF(2)3 with 
derivabive evaluation is not reoommended. Values of the soaled 
algor ith * f various values of <* are found in table 11(b) for the RKT( 2)3 


FOURTH-ORDER SCALED ALGORITHMS FOR THE RKF(3)1* FORMULAS 


. tJO* f f the fourth-order, soaled RK algorithms is quite similar 

RK ^ 2 * 3 methods. Coefficients for a four-stage, fourth- 
al8 °^thm with an embedded third-order method for error estimation 

! re pr ' 8Sented in tabl ® *V(a). The construction of the 
soaied algorithm for these coefficients requires only one additional deriva- 
tive evaluation. Equations (11a), (11b), and (11o) are imposed along with 


8 


1 


• 4 § t §? 

p 5,1 ■ £ 1 ^5,l a i * «5 f 2 


( 16 ) 


* • i 

and Ci ■ 0. The free parameters fer this algorithm are «rj and flg jj. Ail 
ether coefficients are determined from existing eeefflelents using equations 
(8a) or (8b), or from the equations of oonditlon. 

The Cj coefficients may bo wrltton aa 




• 4 • 

Ci+1 0 S (aij / ((J+1)«i + i)) 

JbI 


(17) 


for 1 a 1, 2, 3, and 4, where the aj. j coefficients are determined from 
the Lagrange polynomial: 


MX) 


n* (X -U k ) 

kai 

i; 


a 2 ai j X.1- 1 

J «1 


(18) 


where 

with 


n a 4 and where the defining set M a {p p n ) 
* denoting k it 1* 


. • # * t, 

(a 2l a 3 , a 4 , o 5 ) 


The unknown 85 >k Coefficients may be written as 


65,1 ■ * (i, c * **•> ) ' c 5 


(19) 


« » « • * 

%,2 8 <*2 - K , 03) / Og (<*2 - O3) 


( 20 ) 


and 


9 


( 21 ) 


* « • • • 

05,3 8 ( R g ” K 1 ag) f «3 («3 - ag) 

Tho K*, 1 b 1,2 are defined by 


• • ii » §1 

K i 0 p 5,i " 05,1 «1 - 05,4 <*4 


with 


• • « 

^5,2 8 (1/12 - l Cj Pj f 2 ) f C$ (22) 

J®2 


• j -1 • »2 

P J*2 8 k * 1 0j,k <*k (23) 


Then 05,0 is determined from equation (3). 

• « • 

While the values of 05 and 35^ may be varied to reduce G5 (eq. (9)), 

the system is most sensitive to the 015 parameter. In fact, suoh a small reduo - 

• « 

tion is achieved in G5 by changing 35,4 » that this parameter is set equal to 

zero in all of the scaled RKF(3)4 algorithms. The 05 parameter, however, must 

e 

be seleoted more carefully in order to ensure a nearly minimal value of Gg. 
Values of the soaled coefficients for the RKF(3)4 algorithm, with various 
values of a, are presented in table IV(b). 


FOURTH- AND FIFTH-ORDER SCALED ALGORITHMS FOR THE RK(4)5 FORMULAS 


The complexity of the equations of oondition for the fifth-order, soaled 
algorithm requires that two new derivative evaluations be added. A fourth- 
order, soaled solution may, however, be achieved by the use of only one addi- 
tional function evaluation. Quite often, for intermediate data output, the 
lower-order solution is sufficient and is all that is desired if the defining 
algorithm is being operated as a fourth-order method. (The fourth-order, 
soaled solution for the RK(4)5 will vary slightly from that of the RK(3)4 
formula beaause an additional assumption has been imposed and because five 


10 


derivative evaluations have already been wade for the defining algorithm.) For 

the fourth-order, scaled algorithm, <*6 is the only free parameter, while for 

» » » » 

th© fifth-order, soaled solution Pg g and' Pg § remain as free 

parameters. Coeffiolonta for a six-stage, fifth-order Runge-Kutta algorithm, 
RKF(4)5 II, developed by Fehlberg (ref. 3), are presented in table V(a). 


The simplifying assumptions, imposed during the development of the 
defining algorithm, 


• J-1 « • » 2 

P J,1 8 2 0 J,k a k 8 a J /2 (24) 

k=1 


and 


• j-1 

P J»2 8 2 


• » 2 » 3 

0 J,k a k 8 °j /3 


(25) 


are imposed, along with C y a 0, for J = 2, ...,P« t where P» s 6 for the 
fourth-order solution and P» *7 for the fifth-order solution. 


The Ci ooeffioients may be written as 


• P«-1 « 

c i+1 8 2 ai,j/((j+1) <* i+ i) (26) 

jel 


for i s 1,2, ..., p*-1, 

« « 


on the set M 8 {o^, 03, 


where the ai t j terms are defined by equation (18) 

« 

. . . , } . Then 


Cq 8 



(27) 


For the fourth-order solution, 


11 


( 28 ) 


and 


Then defining 


and 




( 29 ) 


( 30 ) 


( 31 ) 


* • * « « 

06,3 8 < k 3 " K S(«l| ♦ o 5 ) + K-i 04 6(5) /D3 (32) 

where 

„ • • e 1 • e 1 

K i 8 p 6,i - $6,1 <*1 - 06,2 <*2 

with 


• • • • • 

D 3 = 03 (03 - <*4) (03 - 05) 

• § 

where ^ 06,4 an{ * ^06,5 a **® determined by permuting the indices on a 3 , 04, 
and 05. Then 06,0 oan be evaluated from equation ( 3 ). (While equation (29) 
is not requirv for a fourth-order solution, its inolusion helps reduce higher 
order terms used in forming Gg, equation (9).) 
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mmmmm 




For the fifth-order, scaled BKF(4)5 formulas, with the C* coefficients 
determined from equation ( 26 ), 


§ • * § • 

06,1 b (Ki - Ko aj)/C6(a6 - Qt?> 


and 


« • • 

07,1 3 ’ KO ot6)/C 7 


« • 

(a 7 - a6) 


where 


for i s 1 , 2 and 



(33) 


(34) 


(35) 


* * 

where 06 2 is a free parameter. With 06,5 as an additional free parame- 

* # 

ter, the remaining 06 ,k coefficients may be written as 


and 


where 


• * i • * 

06,3 8 <K 2 - K t ot n ) 3 (a 3 -<* 4 ) 


(36) 


« tit* 

06, M 8 < k 2 - K 1 03)^4(014 -a 3 ) 


(37) 
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*i 8 p 6,l - 06,1 <*1 - 06,2«2 " 06, 


Using the definition for Pj , 3 , equation (30), for J s 2,3* 6 


• / 6 « « 1 • 

P 7|3 c l 1/20 - 2 Cj P J(3 J/C 7 
Js2 


ft ft i 

Using the values of P?,i» P 7| 2 and P 7|3 and defining 


• j-1 « • • 

Sj a Cj 2 0,j,k (Jfcj / C 7 

k=3 


for J a 4,5,6, and 7, 


07,6 



2 B J 0j,1 

j»3 


where 


A 3 a 


«• ft ft ft ft ft ft ft 

(K 3 - K2(oui + 05) + Ki 04 ag)/a 3 (a 3 *» ct4>(a 3 - 015) 


and 


b 3 


» • • • • ••»•• • 

-O6(ot6 - 016(04 ♦ 05) ♦ 04 05)/o 3 (n 3 - 04) (o 3 - 04) 


and 


• • »i » »i 

K i = p 7,i “ 07,1 °1 " 07,2 “2 



(B^p U) Then ^ ^ 85 deterroined b y permuting the indioea of the oT'a. 


P7J ■ Aj + Bj 3 7|6 


§ 

f°p J s 3,1) and 5, and 37,0 nay be evaluated from equation (3), 

, The f ' ourth -° rdQr aolutlon la not extremely aenaltlve to the aeleotlon of 

as so that aooeptable error ooefflolenta nay be determined from a given 4 
for a range of 0. For a e (0.00, 0.30], 4 • 0.60, while for o e 

f or ^the°f if th-^der°algori thm ^however |°must M°2ino14 

enaure a near-nlnlnal value of 0* (eq. (9)). The fifth-order solution la most 
sensitive to the eg and a, parameters, and (In general) 4,2 * 4,5 * 0 

X^“rC variations ln\hr^r™urr;4f L^y^o^irslg^anr^anlef 
table V(o)! 088 f ° th fiffch “ order * sealed algorithm are listed in 


EMBEDDED THIRD-ORDER SOLUTIONS AT NO ADDITIONAL DERIVATIVE EXPENSE 

rHS 1 “5?: S"~" at 

S)‘= M S 6 7 0r r ffl0Unt “- ^us 1 whne"h^ a!gorlt^s are :f n0n ' 

ror^r"^ 

f r? 0 ^ . 8 * a low-order estimate of the solution is auffMMonf «.!!?« 

valJe of ®o & a t °no '"add it i onal “der iva t iv^wpen 86?° ** fitenera ^ ed for an * 


15 












A minimum of two nonzero C k coefficients, k > 2, la required to gener- 
ate a third-order, soaled algorithm. The additional derivative evaluations 
used in the RK(3)4 and RK(4)5 formulas allow a ohoioe of one (or two) more 

nonzero Cjj's for oauaing higher-order truncation error ooeffloienta to van- 


ish. However, the norm of the error ooeffiolents G^, Gg, and Gg (eq. (9)) 

£ 

is generally smaller (for the algorithms presented) when only two Cu coeffi- 
cients are nonzero. For the RR(2)3 algorithm 


c i * o(2o - 3aj)/(6 ai((*i - aj)) (41) 

For is2,js3 and is3,Js2. For three nonzero C*, 

« a(3o 2 - 4a(ou ♦ a k ) + 6 oj a k ) 

Ct * ■ (42) 

12 on (a* - aj) (a* - ot k ) 


where (i,j,k) e { (2,3,4), (3,4,2), (4,2,3) ) for the soaled RKF(3)4 algorithm. 
Un general, however, equation (41) should be used for the RKF(3)4 and RKF(4)5 

ft ft ft 

formulas presented, although the nonzero C k 's need not be C2 and C3.) 

A oase of particular interest oocurs for the RKF(4)5 formulas. Using 

* 

equation (42) to generate the C k coefficients as a function of a, the one 
nonzero fourth-order truncation error coefficient may be written as a qua- 
dratic in o. Thus, two values of the parameter a exist that will give a 
fourth-order algorithm. For the RKF(4)5 formulas, a * 0.60 and o = 1.0 give 
a fourth-order, soaled Runge-Kutta algorithm regardless of the defining set of 

§ 

coefficients or the selection of C k coefficients. Because a fourth- and a 
fifth-order solution are already available at o « 1.0, however, the only 
value of interest is a * 0.60. 

« « 

^e c 2 and Co coefficients are determined by equation (41) for the 
hkf( 2)3 algorithm, the fourth-order, scaled truncation error coefficients are 
of the same order of magnitude as those for the defining algorithm, 

0* 8 0.5144E-01. Thus, in using the RKF(2)3 algorithm, no additional deri- 
vative evaluations should be used* Rather, the solution determined by 
equation (41) should be employed. For the RKT(2)3 formula, the fourth-order 


16 







• 

truncation error coefficients are considerably smaller* 04 » Q.1447E-Q3. While 
a solution obtained by using equation (41) Is of third order, a scaled solution 
may not be as accurate as the solution obtained by the defining algorithm that 
has much lower fourth-order error coefficients. While this scaled solution is 

• 

quite accurate near t 0 and t 0 +h (for e a 0,05, G 4 a 0.41588-03 and for 

« • 

0 a 0*95, Gj| a 0.4533E-02), the error terms may be as large as Q|| a 0. 16658- 

01 throughout the interval. However, the scaled solution presented in the sec- 
tion entitled n Thlrd-0rder Scaled Algorithms" is aetually of fourth order. 

Thus, for some problems, the additional derivative evaluation may produoe a 
more acourate scaled solution. The simplloity of the soaled solution offered 
by in equation (41), however, may compensate for the reduoed aoouraoy. 

A study of the RKF(3)4 algorithm shows that equation (41), defined for 

• • 

Cg and Cj| , produces the most aeourate thied-order, soaled solution for 

« • « 

o.< 0.65. However, for o >0.65, using a Cg, C 3 solution with C 4 e 0 , 

# 

will give the best results. (Eq. (42) produoes a G 4 that is approximately 
twice the magnitude of the best solution obtained when using eq. (41).) 

• • 

For 0 < 0.05, O 4 _< 0.3183E-03. Throughout the interval, however, 04 may be 
as large as 0.9390E-02. 

The most important feature of the low-order, soaled RKF(4)5 formulas is 
the existence of a fourth-order solution at 0 a 0.60, whioh requires no addi- 

« 

tional derivative- evaluation. This fourth-order solution with the Cu's, de- 

§ • • * 

fined by equation (42) (using Cg, C 3 , and C 5 )-, gives 0 5 a 0.1978E-02. For 
other values of o, a good third-order, soaled solution may be obtained using 

« « 

equation (41) for Cg and Co for all values of a, although for o e (0.35, 
0.55), a slightly better solution may be obtained using equation (41) to 
* * » • 
determine Cg and C 5 . For a < 0.05, G 4 < 0.2884E-03 with Gjj < 0.5376E-02 

§ “ 
for o < 0.75. For o > 0.75, G 4 < 0.2080E-01* 


APPLICATION OF THE SCALED ALGORITHM 


The most direct manner of applying the soaled Runge-Kutta algorithm is to 
demand data output at a speoifio value of the independent variable. With some 
maneuvering, however, numerioal integration oodes oan be written that will 
output information by specifying values of a dependent variable or by satis- 
fying a constraint equation. While the scaled algorithms have been developed 


for a wide range of a {a a Q.Q5K, K * T, ..., 19), data will generally be 
requested for values of the Independent variable t s t Q + oh, for values of 
a that wore not studied. The simplest manner of applying tho algorithm would 
be to reduce the original stop size slightly so that the requested output point 
corresponds to a t « t 0 + uh for a o of a determined algorithm. For exam- 
ple, if tho solution is required at t ■ t 0 ♦ 0.87 hi (hi being the predicted 
step size), tho user would let the step size bo h b 0.87 hi /0 . 90, which la 
less than a 4-peroent reduction from the predlotod step size. For output 
near t 0 , the percentage of step-size reduction may bo as large as 50 peroent 
with no soaled output possible for a < 0.05 without generating the ooeffi- 
oients from the sealed algorithm formulas. These problems near t s t 0 may be 
skirted in a number of ways, for example, by interpolation or by using the 
third-order algorithm that requires no derivative evaluations, whioh is valid 
for any ohoioe of a and whioh has low fourth-order trunoatlon error coeffi- 
cients for or near zero. 

Difficulties in applying the algorithm in suoh a manner may arise (1) if 
the value of the independent variable for data output is not known before the 
step is taken, (2) if storage limitations restriot the number of soaled coeffi- 
cients that may be retained, or (3) if dense output is requested within the in- 
tegration step. In suoh oases, an interpolating polynomial may be generated, 
using points determined by evaluating the soaled algorithm for several values 
of o. Sinoe the solution and its derivative are known at t 0 and t 0 ♦ h, a 
fourth-order polynomial may be determined by applying a soaled algorithm at 
only one point. For the RKF(4)5 algorithm, a fourth-order solution is avail- 
able for a a 0.60 requiring no additional derivative evaluations. Thus, an 
interpolating polynomial may be generated quite inexpensively. (If output is 
dense in one region of the interval, a different interpolating polynomial gen- 
erated by scaled solutions in the neighborhood of this region would be more ac- 
curate. However, these points oan be generated inexpensively, oompared to re- 
peated Runge-Kutta solutions at small step sizes. ) The disadvantage of using 
the interpolating polynomial is that the solution is subjeot to errors from the 
polynomial as well as from the Runge-Kutta solution; whereas, the scaled solu- 
tion is equivalent to a Taylor series solution of corresponding order. 

While the soaled algorithm is designed to determine information at spec- 
ified values of the independent variable, the methods may be adapted to atop 
the integration by specifying values of a dependent variable or by satisfying 
a constraint equation. The usual technique for stopping the integration on a 
specified value of a function is to bracket the desired solution and to iterate 
until the value of the function is reached. In order to be operated effi- 
ciently, however, a Runge-Kutta method must take large step sizes (consistent 
with truncation error oriteria). While using a Runge-Kutta method to trap a 
solution with smaller and smaller step sizes will produce an aoourate solution, 
the expense of computing the derivative evaluations may be prohibitive. The 
use of soaled algorithms to insert data points allows the formulation of an 
Interpolating polynomial. This polynomial may be analyzed efficiently to 
determine the value of the independent variable* whioh causes the constraining 
funotion to be satisfied. If the accuraoy of the Runge-Kutta solution is 
required; i.e., if the inaccuracies introduced by the interpolating polynomial 
produce an unacceptable solution, an estimate may be used to determine the 
next step size in an attempt to satisfy the constraining funotion. The 
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TABLE I.« COEFFICIENTS POR RKF(2)3 ACCORDING 
TO FEHLBERO (REP. 3) 



EUCLIDEAN NORM OF TRUNCATION ERROR COEFFICIENTS (TEC), Oil s 0.5144159D-01 

STABILITY LIMITS R: (-0.245M4D+01 , O.O) 

I: (0.0, 0.17908D+01) 


TABLE II. « COEFFICIENTS FOR RKT(P)3 
(a) Unsealed algorithm aooordlng to Betti a 1 



EUCLIDEAN NORM OF TEC O 4 a 0 . 1 M46759D-03 

STABILITY LIMITS R: (-0.2791000D+01 , 0) 

I: (0, 0.2827170D+01) 


Private communication. 



- ■ ■,,< 



0 ^ , o 








TABU! II,-. Continued 
(b) 3oaled algorithm 


i 

| 


4 

SIQMA a 0 , SO0OOOD-O 1 j A (A) b 3,/S. 


CO 8 
Cl a 
C8 8 
C3 8 
C4 a 
B40 a 
B41 a 
B42 8 
B43 s 


419 / g 4QQ 


37 / 56 400 
- 17 / 232 800 

45 125 / 54 708 

271 891 / 475 000 

- 74 533 / 12 996 000 

13 948 237 / 324 900 000 

- 4 559 / 475 000 


SCALED NORM OP TEC 
SCALED STABILITY LIMITS 


Q*(5) a 0.71 562 14D-05 
R« (-0.9478000D+01, 0) 
I* (0, 0, 1 175900D+02) 


SIGMA a 0. 100000D+00; A(4) a 3./5. 


CO a 
Cl a 
C2 a 
C3 a 
C4 a 
B40 a 
B41 a 
B42 a 
B43 a 


109 / 600 
0 

17 / 6 600 

-7/28 200 
3 375 / 4 136 

30 691 / 56 250 
- 1 591 / 145 800 
525 007 / 18 225 000 
- 517 / 28 125 


SCALED NORM OP TEC 
SCALED STABILITY LIMITS 


G*(5) a 0.5347556D-04 
R: (-0.6646000D+01 , 0) 
(0, 0.8300200D+01) 


1 
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TABUS II«» Continued 
(b) Continued 


SXQMA s C, 150000D+00; A(4) » 3 


i 

CO 8 

451 / 2 400 

m, 

Cl 8 

0 

§ 

C2 8 

93 / 16 400 


C3 s 

- 33 / 72 800 


C4 8 

36 125 / 44 772 

t - 

B40 a 

220 857 / 425 000 

Pi- 

B41 a 

- 53 521 / 3 468 000 

B 

B42 a 

10 586 569 / 86 700 000 


B43 a 

- 11 193 / 425 000 

&■ 

SCALED NORM OP TEC G*(5) a 0. 1683013D-03 

B 

SCALED STABILITY 

LIMITS R: (-0.5412667D+01 , 0) 

& 


I: (0, 0.6766133D+01) 


0.200000D+00 : A(4) 


CO 8 
Cl s 
C2 8 
C3 8 
C4 8 
B40 s 
B41 s 
B42 s 
B43 s 


29 / 150 
0 

14 / 1 
- 1 / 1 
500 / 627 
091 / 6 


1 425 
1 650 


3 091 / 6 250 

- 2 761 / 144 000 

568 993 / 3 600 000 

- 209 / 6 250 


SCALED NORM OP TEC Q*(5) a 0.37 1388 ID-03 

SCALED STABILITY LIMITS R: (-0.4690500D+01 , 0) 

I: (0, 0.5851 150D+01) 






TABLE II.- Continued 


(b) Continued 


SIGMA b 0.250000D+00; A(4) a 3./S. 

CO 8 

19 / 96 

Cl B 

0 

C2 b 

5 t 336 

C3 b 

- 1 / 1 632 

C4 s 

375 / 476 

B40 b 

1 411 / 3 000 

B41 b 

- 7 121 / 324 000 

B42 b 

12 397 / 64 800 

B43 b 

- 119 / 3 000 

SCALED NORM OP TEC 

G*(5) b 0.6741268D-03 

SCALED STABILITY LIMITS R: (-0.4207200D+01 , 0) 


It (0, 0.5226880D+01) 

SIGMA s 0.300000D+00; A(4) b 3-/5. 

CO s 

121 / 600 

Cl 8 

0 

C2 s 

33 / 1 600 

C3 = 

- 3 / 8 200 

C4 s 

6 125 / 7 872 

B40 8 

19 557 / 43 750 

B41 s 

- 869 / 36 750 

B42 8 

101 803 / 459 375 

B43 * 

- 984 / 21 875 

SCALED NORM OF TEC 

G*(5) b 0.1080755D-02 

SCALED STABILITY LIMITS R: (-O.3058333D+O1, 0) 


It (0, 0.4766600D+01) 


25 


TABLE II.- Continued 
(t>) Continued 


SIGMA 

b 0.350000D+00; A(4) b 3./5. 

CO s 

491 / 2 400 

Cl B 

0 

C2 b 

931 / 34 800 

C3 *> 

49 / 189 600 

C4 s 

21 125 / 27 492 

B40 s 

138 013 / 325 000 

B4l a 

- 29 099 / 1 216 800 

B42 b 

37 812 607 / 152 100 000 

B43 a 

- 16 037 / 325 000 

SCALED NORM OP TEC G»(5) = 0.1589552D-02 

SCALED STABILITY LIMITS Rj (-0.3594286D+01 , 0) 


I: (0, 0.4409743D+01) 

SIGMA 

8 0.400000D+00; A(4) s 3-/5. 

CO 8 

31 / 150 

Cl 8 

0 

C2 s 

32 / 975 

C3 8 

2 / 1 425 

C4 s 

375 / 494 

B40 8 

3 781 / 9 375 

B41 s 

- 1 817 / 81 000 

B42 8 

550 433 / 2 025 000 

B43 8 

- 494 / 9 375 

SCALED NORM OP TEC G*(5) s 0.219401 ID-02 

SCALED STABILITY LIMITS R: (-O.330725OD+O1 , 0) 


I: (0, 0.4123225D+01) 
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TABLE II.- Continued 
(b) Continued 


SIGMA s 0 . 450000D+00 ; A(4) b 3./5. 


CO B 
Cl 8 
C2 b 
C3 B 
C4 s 
B40 s 
B41 b 
B42 b 
B43 b 


499 t 2 400 

0 

351 / 9 200 

189 / 58 400 

15 125 / 20 148 

105 339 / 275 000 

- 9 061 / 484 000 

3 516 493 / 12 100 000 

- 15 111 / 275 000 


SCALED NORM OF TEC G*(5) s 0.2883997D-02 

SCALED STABILITY LIMITS R: (-0.3222000D+01, 0) 

I: (0, 0.38873330+01) 


SIGMA a 0.500000D+00; A(4) b 1./2. 


CO s 

1 / 

6 

Cl B 

0 


C2 b 

1 / 

6 

C3 B 

0 


C4 S 

2 / 

3 

B40 b 

5 / 

16 

B41 8 

- 1 / 

8 

B42 s 

3 / 

8 

B43 s 

- 1 / 

16 


SCALED NORM OF TEC G*(5) b 0.3654457D-02 

SCALED STABILITY LIMITS R: (-0.3088000D+01 , 0) 

I: (0, 0.3689480D+01) 


TABLE II.- Continued 
(b) Continued 


SIGMA b 0.550000D+00; A(4) a 1./2. 


CO 8 
Cl 8 
C2 a 
C3 8 
C4 8 
B40 s 
B 41 s 
B42 s 
B43 n 


209 / 1 200 

0 

121 / 600 
121 / 34 800 
18 / 29 
1 421 / 4 800 

- 4 279 / 31 104 

317 251 / 777 600 

- 319 / 4 800 


SCALED NORM OP TEC 
SCALED STABILITY LIMITS 


G*(5) s 0.444375 ID-02 
R: (-0.2978545D+01, 0) 
IJ (0, 0.3521364D+01) 


SIGMA s 0.600000D+00; A(4) s 1./2, 


CO 8 
Cl 8 
C2 s 

C3 * 
C4 8 
B40 s 
B41 s 
B42 s 
B43 s 


9/50 

0 

6/25 

3 / 350 

4 / 7 

7/25 

- 173 / 1 152 

12 677 / 28 800 

-7/100 


SCALED NORM OP TEC 
SCALED STABILITY LIMITS 


G*(5) s 0. 5264721 D-02 
R: (-0.2889667D+01, 0) 
!•* (0, 0.3377233D+01) 


TABLE II.- Continued 
(b) Continued 


SIGMA s 0.650000D+00} A(4) a 1./2. 

CO 8 

221 / 1 200 

Cl B 

0 

C2 b 

169 / 600 

C3 b 

169 / 10 800 

C4 8 

14 f 27 

B40 b 

423 / 1 600 

B41 8 

- 1 313 / 8 064 

B42 s 

95 069 / 201 600 

B43 8 

- 117 / 1 600 

SCALED NORM OP TEC 

G*(5) s 0.61051 16D-02 

SCALED STABILITY LIMITS R: (-0.2818000D+01, 0) 

I: (0, 0.3253138D+01) 


SIGMA s 0.700000D+00; A(4) = 1./2. 


CO 8 

14 / 75 

Cl 8 

0 

C2 s 

49 / 150 

C3 s 

49 / 1 950 

C4 s 

6/13 

B40 s 

299 / 1 200 

B41 8 

- 455 / 2 592 

B42 s 

32 543 / 64 800 

B43 * 

- 91 / 1 200 

SCALED NORM OF TEC 

G*(5) s 0.6957139D-02 

SCALED STABILITY LIMITS R: (-0.2761857D+01 , 0) 

Ir (0, 0.3146329D+01) 



CO s 
Cl s 
C2 = 
C3 s 
C4 s 
B40 s 
B4l s 
B 42 s 
B43 8 


TABLE II*- Continued 
(b) Continued 


0.750000D+00; A(4) = 3./10 


51 / 88 

15 / 496 

250 / 1 023 

837 / 8 000 

- 30 079 / 48 000 

9 119 / 9 600 

- 1 023 / 8 000 


SCALED NORM OP TEC G* ( 5 ) s 0.7719994D-02 

SCALED STABILITY LIMITS R: (-0.2720267D+01 , 0) 

I: (0, 0.3055080D+01) 

SIGMA s 0.800000D+00: A(4) s 3./10. 


CO 8 
Cl 8 
C2 s 
C3 8 
C4 8 
B40 s 
B41 s 
B42 s 
B43 s 


608 / 975 


125 / 741 


1 425 


43? / 
- 14 473 / 


6 250 
18 000 


536 497 / 450 000 
- 494 / 3 125 


oSJHS! 0F TEC Q*(5) * 0.840 109 ID-02 
SCALED STABILITY LIMITS R: (-0.2693250D+01 , 0) 

I: (0, 0.2978425D+01) 




TABLE II— Continued 
(b) Continued 


SIOMA ■ 

0.6S0000D+00) A(4) a 3. /IQ. 

CO e 

203 / 1 200 

Cl B 

0 

C2 b 

19 363 / 29 400 

C3 B 

12 427 / 178 800 

C4 b 

750 / 7 301 

B40 s 

9 983 / 600 000 

B41 b 

- 283 883 / 259 200 

B42 b 

51 368 611 / 32 400 000 

B43 b 

- 124 117 / 600 000 

SCALED NORM OP TEC G»(5) s 0. 996999 3D-02 

SCALED STABILITY LIMITS R: (-0.2681647D+01 , 0) 

I: (0, 0.2916259D+01) 

SIOMA 

b 0.900000D+00; A(4) b 1./10. 

CO 8 

7/50 

Cl 8 

0 

C2 8 

1 431 / 2 050 

C3 s 

216 / 2 275 

C4 s 

250 / 3 731 

B40 s 

- 6 643 / 50 000 

B41 s 

- 3 991 / 3 000 

B42 s 

134 029 / 75 000 

B43 = 

- 11 193 / 50 000 


SCALED NORM OP TEC Q*(5) s 0.9706279D-02 

SCALED STABILITY LIMITS R: (-0.2688444D+01, 0) 

I: (0, 0.286931 1D+01) 


TABUS II.- Concluded 
(b) Concluded 


SIGMA b 0.9500000+00} A (4) 



CO 8 

197 t 1 200 

Cl 8 

0 

C2 b 

11 191 / 16 200 

C3 8 

27 797 / 217 200 

C4 s 

250 / 14 661 

B40 a 

- 74 753 / 200 000 

B4t a 

- 45 049 / 16 000 

B42 a 1 

501 437 / 400 000 

B43 « 

- 92 853 / 200 000 

SCALED NORM OF TEC 

G*(5) 8 0.1094176D-01 

SCALED STABILITY LIMITS R: (-0.2718947D+01 , 0) 


Is (0, 0.2839242D+01) 
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TABLE III. - Concluded 


1 / P g 1-1 1 

*.« ‘ i [il Cl >i Bl)J “ J ' 5 

T5 ' 3 ’j(il° l 01 J.1 6l,J “ ?J ' » 


t5 ' 4 'il°‘ 0,1 >a 6lJ kli Bj,k ° k ’ 55 


1 / P i*1 Q 1 

* 5 (i Cl jJ, Bl ' J ’ S 


T5 - 6 ' ll Cl \i 6iJ “ J to, Sj,k ° li 5o 


, / p i-1 J-1 2 1 

T5 ' 7 ’i(u 3 Cl A 6l ' J to. 


T5 - 8 \*4 Cl M BlJ krf Bj,k U 1 BM ° l ‘^0 

l5 ' 5 ■ 2 ( il Cl ‘!f, Bl,J OJ>2 ‘ io ) 
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TABUS IV. - COEFFICIENTS FOR RKF(3)*4 
(a) Urmoaled algorithm 



EUCLIDEAN NORM OF TEC 
STABILITY LIMITS 


0(5) n 0.7141547D-02 

R: (-0.4109200D+01, 0) 
It (0, 0.3214450D+01) 
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TABkE IV. - COEFFICIENTS FOR SCALED RKF(3)4 
(b) Scaled algorithm 


SIGMA s 0.8QGOG0D-01} A(8) a p,/3. 


CO i 
Cl 8 
C2 s 
C3 a 
CD e 
C5 b 
B50 c 
B51 b 
85a = 
B53 b 
BSD b 


2 072 939 / 11 536 000 
0 

- 12 393 / 03 200 000 

1 495 823 / 12 954 240 000 

- 3 911 / 69 600 000 

150 651 / 200 680 

29 852 683 / 47 574 000 

1 798 / 330 375 

1 063 221 / 22 906 000 

- 3 924 851 / 309 231 000 

0 


SCALED NORM OF TEC 
SCALED STABILITY LIMITS 


G*(5) s 0.5603873D-05 
8? (-0.10114000*02, 0) 
I: (0, 0.12569000*02) 


SIGMA 8 0. 100000D*00$ A(5) a 16. /25. 


CO 8 
Cl s 
C2 a 
C3 * 
C4 b 
C 5 s 
B50 a 
B51 a 
B52 a 
B53 8 
B54 b 


1 152 549 
- 333 739 
167 

20 663 828 
4 463 312 576 
- 27 858 100 608 
3 330 510 639 192 
2 456 985 733 216 


517 603 / 2 304 
0 

/ 890 240 
/ 1 299 
/ 1 755 
/ 26 690 
/ 7 748 
/ 2 582 
/ 33 578 
t 100 736 
0 


125 

104 


000 

000 

168 000 
000 

107 392 
935 546 875 
978 515 625 
720 703 125 
162 109 375 


Scaled norm op tec 

SCALED STABILITY LIMITS 


0*(5) e 0.4156062D-04 
R: (-0.70950000*01, 0) 

I* (0, 0.88678000*01) 


TABLE 1V*« Continued 
(b) Continued 

SIGMA m 0. 150QQQD±QQ ; A(5) a 16. /2G 


GO a 



43 

261 

/ 

192 

000 

Cl 8 





0 

C2 8 


2 

184 

813 

/ 

1 

630 720 000 

C3 a 



362 

551 

/ 

554 

112 000 

C4 b 

5 

838 

- 6 

687 

/ 

18 

080 000 

C5 ■ 

359 

375 

t 

7 

552 290 816 

B50 a 

393 

604 

655 

104 

/ 

729 

794 921 875 

B51 s 

1 

855 

138 

176 

/ 

729 

794 921 875 

B52 b 

1 241 

843 

094 

576 

/ 

9 

487 333 984 375 

B53 s 
B54 8 

- 310 

926 

657 

216 

t 

0 

9 

487 333 984 375 

SCALED 1 

NORM OP 

TEC 


G*(5) s 

0.1296665D-O3 

Scaled i 

STABILITY LIMITS 

1 Rt 

( 

-0.5783333D+01. 0) 





It 

(0, c 

I.7225667D+01) 


0.200000D+00} A(5) a 1./8. 


Cu s 
Cl 8 
C2 a 
C3 8 
C4 a 
C5 a 
B50 s 
B51 a 
B52 a 
B53 a 
B54 s 


2.9 671 

- 17 356 

42 

48 249 
16 014 162 

- 56 240 
73 783 218 
15 586 385 


897 / 180 
0 

029 / 196 
829 / 272 
103 / 1 

856 / 20 
551 / 72 
379 / 188 
033 / 313 
647 / 470 
0 


300 000 
610 000 
462 500 
582 055 
374 784 000 
476 000 
624 064 000 
436 096 000 


SCALED NORM OP TEC 0*(5) 8 0.7441866D-03 

SCALED STABILITY LIMITS Rj (-0.5017500D+01, 0) 

It (0, 0.6246150D+01 ) 





TABLE IV.- Continued 
(b) Continued 


3IUMA a 0.250QCQB+QQJ A(5) a 3./5. 


CO B 
Cl B 
C2 B 

C3 s 
C4 b 
C5 « 
B50 b 
B51 « 
B52 b 
B53 * 
B54 * 


11 057 / 55 296 
0 

124 659 / 7 055 360 

2 401 / 11 860 992 
-9/21 760 
628 250 / 802 791 
3 372 873 / 7 180 000 

- 431 919 / 12 565 000 

277 074 513 / 1 306 760 000 

- 8 851 689 / 186 680 000 

0 


SCALED NORM Of TEC 
SCALED STABILITY LIMITS 


Q»(5) s 0.5096777D-03 
R: (-0.4507200D+01, 0) 
I; (0, 0.5578000D+01) 


SIGMA = 0.300000D+00; A(5) ■ 3»/5* 


CO 8 
Cl 8 
C2 b 
C 3 = 
C4 s 
C5 s 
B50 b 
B 51 s 
B52 a 
B53 * 
B54 a 


57 851 / 288 000 
0 

4 310 577 / 247 520 000 
304 927 / 49 296 000 
- 2 637 / 820 000 

$ 402 875 / 6 937 938 

9 318 525 453 / 21 611 500 000 

2 632 446 / 675 359 375 

60 133 894 317 / 280 949 500 000 

- 6 900 061 371 / 140 474 750 000 

0 


SCALED NORM OP TEC 
SCALED STABILITY LIMITS 


G*(5) s 0. 808207 4D-0 3 
R: (-0.4140333D+01, 0) 
I: (0, 0.5085633D+01) 


TABLE IV.- Continued 
(b) Continued 


SIOMA a 0.350000D+Q0} A(5) b 9./16, 


CO a 

19 065 713 / 

103 680 000 

Cl a 

0 


C2 a 

3 252 504 213 / 

59 321 600 000 

C3 a 

28 118 111 / 

6 514 560 000 

C4 a 

- 536 403 i 

205 600 000 

C5 a 

1 634 729 984 / 

2 152 165 545 

B50 a 

5 403 228 165 44 1 / 

13 077 839 872 000 

B51 a 

- 49 882 382 361 / 

817 364 992 000 

B52 a 

45 186 311 217 159 / 

170 Oil 918 336 000 

B53 a 

- 2 355 259 443 201 / 

42 502 979 584 000 

B54 a 

0 


SCALED NORM OP TEC G*(5) a 0. 1 172146D-02 

SCALED STABILITY LIMITS Hi (-0.3864000D+01., 0) 


I: (0, 0.4704371D+01) 


SIGMA a 0.400000D+00; A<5) a 11./20. 


CO a 
Cl a 
C2 a 
C3 a 
C4 a 
C5 a 
B50 a 
B51 a 
B52 a 
B53 a 
B54 a 




446 

/ 

A 

2 

in 

=* 




24 

057 

u 

/ 

328 

250 




16 

807 

/ 

1 

304 

550 



- 94 

/ 

14 

625 



21 

452 

000 

/ 

28 

987 

101 


1 226 

671 

061 

/ 

3 

217 

800 

000 

- 5 

764 

473 

/ 

134 

075 

000 


372 

592 

737 

/ 

1 

394 

380 

000 

- 2 318 

720 

327 

/ 

41 

831 

400 

000 


0 


SCALED NORM OP TEC G*(5) a 0. 1594736D-02 

SCALED STABILITY LIMITS R: (-0.3649500D+01, 0) 

It (0, 0.4398850D+01) 


TABLE IV.- Continued 
(b) Continued 


SIGMA 

a 0.450000D+00; A(5) a 1./2. 

CO a 

208 927 / 

1 280 000 

Cl a 

0 


C2 8 

1 010 269 341 / 

6 572 800 000 

C3 8 

6 043 317 / 

490 880 000 

C4 a 

- 4 941 / 

800 000 

C5 S 

15 ?76 / 

23 305 

B50 a 

48 018 463 / 

126 208 000 

B51 a 

- 1 170 873 / 

7 888 000 

B52 a 

541 971 837 / 

1 640 704 000 

B53 8 

- 1 598 723 / 

25 636 000 

B54 8 

0 



SCALED NORM OP TEC G*(5) a 0.2060460D-02 

SCALED STABILITY LIMITS R: (-0.3480222D+01, 0) 



I: (0, 0.4148133D+01) 


SIGMA a 0.500000D+00; A(5) a 12. /25. 

CO a 

2 201 / 13 824 

Cl 8 

0 

C2 a 

6 561 / 33 280 

C3 a 

45 619 / 1 617 408 

C4 8 

-1/80 

C5 8 

78 125 / 124 416 

B50 3 

138 009 / 390 625 

B51 a 

- 52 704 / 390 625 

B52 a 

1 635 471 / 5 078 125 

B53 a 

- 306 936 / 5 078 125 

B54 a 

0 

SCALED NORM OF TEC G*(5) a 0.2556559D-02 

SCALED STABILITY LIMITS R: (-0.3345400D+01, 0) 

Is ( 0 , 0.3938940D+01) 




TABLE IV.- Continued 
(b) Continued 


SIGMA ■ 0.55Q0Q0D+00} A(5) b 2./5. 


CO 8 
Cl B 
C2 b 
C 3 s 
C4 b 
C 5 a 
B50 s 
B51 s 
B52 s 
B53 b 
B54 8 


313 009 

563 920 137 
99 648 703 
- 226 391 
5 384 875 

6 165 113 549 
- 227 532 591 

7 504 919 433 
- 15 211 290 941 


2 304 000 


1 680 
3 339 
18 720 
10 540 

16 154 
673 109 

17 500 
210 010 


640 000 
648 000 
000 
764 

625 000 
375 

843 750 
125 000 


SCALED NORM OP TEC G«(5) s 0. 306221 4D-02 

SCALED STABILITY LIMITS R: (-0.323, <4D+01, 0) 

Is (0, 0,376^821*01) 


SIGMA s 0.600000D+00; A(5) s 1./3, 


CO 8 

293 / 

2 500 

Cl 8 

0 


C2 8 

1 554 957 / 

3 575 000 

C3 = 

69 629 / 

1 495 000 

C4 s 

- 843 / 

50 000 

C5 s 

8 463 / 

20 240 

B50 s 

174 374 / 

453 375 

B51 8 

- 159 072 / 352 625 

B52 s 

8 735 067 / 

18 336 500 

B53 s 

- 1 804 649 / 

23 575 500 

B54 s 

0 



SCALED NORM OP TEC 
SCALED STABILITY LIMITS 


G*(5) * 0.3559569D-02 
Rs (-0. 3155000D+01, 0) 
Is (0, 0. 3613317D+01 ) 


TABLE IV.- Continued 


(b) Continued 



SIQMA b 0.650000D+00; 

A(5) b 1 

,/4. 


CO b 


111 

507 

/ 

1 280 

000 


Cl a 




0 




C2 8 

677 

766 

609 

t 

1 299 

200 

000 

C3 = 

52 

094 

497 

/ 

746 880 

000 


C4 b 

- 3 

624 

881 

/ 

160 800 

000 


C5 b 

27 

300 

352 

/ 

79 361 

835 


B50 s 

240 021 

016 

229 

/ 

655 208 

448 

000 

B51 * 

- 7 210 

754 

803 

/ 

13 650 

176 

000 

B52 8 

106 401 

934 

989 

/ 

218 402 

816 

000 

B53 = 

- 12 231 

119 

663 

/ 

163 802 

112 

000 

B54 s 

4 

l 


0 




SCALED NORM OF TEC 


G*(5) a 0.4032702D-02 

SCALED STABILITY LIMITS R: 

: (-0.3092769D+01 , 0) 




l! 

» (0, 0. 34874 15D+01) 


SIGMA s 0.700000D+00] 

! A(5) a 1 

I./12. 

CO 8 

_ 

196 

439 

/ 

1 440 

O 

O 

O 


Cl a 




0 




C2 a 

880 

415 

487 

/ 

1 445 

600 

000 

C3 * 

554 

244 

439 

/ 

6 280 

560 

000 

C4 a 

mm 

792 

673 

/ 

33 900 

000 


C5 s 

24 

446 

592 

/ 

52 696 

985 


B50 8 

15 469 

714 

361 

/ 

73 339 

776 

000 

B51 s 

- 406 

078 

043 

/ 

1 018 

608 

000 

B52 s 

43 

531 

047 

/ 

137 936 

500 


B53 a 

CM 

in 

CM 

1 

411 

309 

/ 

953 417 

088 

000 

B54 a 




0 




SCALED NORM OF TEC 


Q* (5) b 0.4462969D-02 

SCALED STABILITY LIMITS R 

• 

« 

(-0.3051000D+01, 0) 




I 

• 

♦ 

(0, 0.3382686D+01) 


TABLE IV. - Continued 
(b) Continued 


SIOMA 

a 0.750000D+00; A(5) a 1./25. 

CO 8 

- 2 453 / 6 144 

Cl 8 

0 

C2 s 

30 292 137 / 49 653 760 

C3 8 

1 063 643 / 7 707 648 

C4 8 

- 4 347 / 124 160 

C5 s 

14 375 000 / 20 948 799 

B50 s 

379 127 171 / 4 600 000 000 

B51 8 

- 37 545 111 / 287 500 000 

B52 8 

6 133 430 349 / 59 800 000 000 

B54 8 

0 


SCALED NORM OF TEC G«(5) s 0.4985026D-02 

SCALED STABILITY LIMITS R: (-0.3030400D+01, 0) 



I: (0, 0.32984130+01) 


TABLE IV.- Continued 
(b) Continued 

SIGMA 

s 0.800000D+00; A(5) a 1./25. 

CO 8 

- 79 / 250 

Cl 8 

0 

C2 s 

553 311 / 942 500 

C3 8 

340 942 / 1 759 875 

C4 s 

- 10 336 / 226 875 

C5 s 

8 828 125 / 15 200 988 

B50 8 

96 032 954 / 3 310 546 875 

B51 8 

28 671 192 / 1 103 515 625 

B52 8 

- 243 083 133 / 14 345 703 125 

B53 * 

84 128 884 / 43 037 109 375 

B54 s 

0 

SCALED NORM OF TEC G*(5) 8 0.6022232D-02 

SCALED STABILITY LIMITS R: (-0.3034500D+01 , 0) 

I: (0, 0. 32352 12D+01) 


TABLE IV. - Continued 
(b) Continued 



SIOMA 8 0.850000D+00} A(5) 

S 

1./25. 

CO 


-603 

317 

/ 

2 

304 

000 


Cl * 




0 




C2 s 

17 122 

676 

913 

/ 

30 

734 

080 

000 

C3 s 

10 521 

438 

907 

/ 

43 

145 

856 

000 

C4 s 

- 11 

352 

209 

/ 

231 

640 

000 


C5 * 

3 931 

093 

750 

/ 

7 

710 

428 

943 

B50 a 

- 722 270 

252 

129 

/ 

23 

586 

562 

500 000 

B51 s 

394 747 

987 

977 

/ 

1 

965 

546 

875 000 

B52 8 

- 30 750 068 

734 

659 

/ 

204 

416 

875 

000 000 

B53 8 

12 397 885 

510 

507 

/ 

613 

250 

625 

000 000 

B54 a 




0 




SCALED NORM OF TEC 


0*(5) 8 

0.7423369D-02 

SCALED STABILITY LIMITS R: 

: ( 

-0.3071 176D+01, 0) 




I! 

: (0, 0.3195200D+01) 


SIGMA s 0.900000D+00; 

! A(5) 


I./25. 

CO 8 


- 7 

279 

/ 

32 

000 



Cl 8 




0 





C2 8 

1 000 

703 

403 

/ 

1 

911 

520 

000 

C3 8 

28 

214 

151 

/ 

99 

632 

000 


C4 s 

- 

202 

743 

/ 

4 

820 

000 


C5 8 

98 

203 

125 

/ 

212 

176 

882 


B50 s 

- 33 318 

422 

659 

/ 

392 

812 

500 

ooo 

B51 s 

1 471 

910 

874 

/ 

4 

091 

796 

875 

B52 s 

- 462 500 

813 

457 

/ 

1 

702 

187 

500 000 

B53 s 

93 979 

832 

093 

/ 

2 

553 

281 

250 000 

B54 8 




0 



% 

SCALED NORM OF TEC 


G*(5) 8 

0.8875535D-02 

SCALED STABILITY LIMITS 

1 R: 

( 

-0.3159222D+01, 0) 




I: 

(0, c 

I.3181233D+01) 
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TABLE XV.- Concluded 
(b) Concluded 


SIGMA s 0.9500000+00; A(5) s 1./10, 


CO 8 
Cl B 
C2 b 
C3 b 
C 4 s 
C5 s 

650 b 
B51 s 
B52 s 
B53 = 
B54 s 


3 273 / 256 000 
0 

378 173 853 / 805 120 000 
129 147 389 / 409 728 000 
- 2 222 677 / 86 880 000 
83 048 000 / 364 430 649 
- 70 545 13? 063 ✓ 284 736 000 000 

40 552 526 367 / 41 524 000 000 

6 270 301 089 333 / 8 636 992 000 000 

89 886 853 133 / 925 392 000 000 

0 


SCALED NORM OP TEC 
SCALED STABILITY LIMITS 


0*( 5) a 0.1004420D-01 
R: (-0.3354842D+01, 0) 
I: (0, 0.3192989D+01) 


TABLE V.- COEFFICIENTS FOR RKF<4)5 
(a) Unaoaled aix-atage, fifth-order algorithm 


! 



EUCLIDEAN NORM OF TEC 


06 * 0.3355745D-02 


STABILITY LIMITS 


Rs (-0.3677700D+01, 0) 

I: (0.2046040D+01 , 0.3606560D+01) 


4 
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TABLE V Continued 
(b) Soalo 1» fourth-order algorithm 


SIGMA 


a 0.05DQ} A( 6 ) s 0.60D0 


C( 0 ) 

C( 1 ) 

C( 2 ) 

C(3) 

C(4) 

C(5) 

C( 6 ) 

B(6,0) 

B( 6 , 1 ) 

B( 6 , 2 ) 

B( 6 , 3 ) 

B(6,4) 

B(6,5) 

SCALED 

SCALED 


8 23615711/129600000 

s 0 

8 2738 14V553078125 
8 6812112671/10482436800000 
a -218411/582000000 
» -481016/155100000 
8 5362046875/6573904758 
8 1 49 1 ^9243/857927500000000 

8 5178062961/53620468750 
8 -29286 158641 302 / 3 1837 153320 “? 1 25 

: 1f?L 2133463595031/ ^^ 

- 95128259521667/1072409375000000 

NObK ?S^ 6M27/5898S5 ’56250000 
«OBM_°P TEC, 0*(5) = 0.1B27626D-0B 

(-0.6246000+01, 0 ) 
I} ( 0 , 0.3606560D+01) 


STABILITY LIMITS 


SIGMA 8 0.10D0; A( 6 ) s 0.60D0 


C( 0 ) 

C( 1 ) 

C( 2 ) 

C(3) 

C(4) 

C(5) 

C( 6 ) 

B( 6 , 0 ) 

8 ( 6 , 1 ) 

B(6,2) 

B(6,3) 

B(6,4) 

8(6,5) 

SCALED 

SCALED 


155536 1/8100000 
0 


8 7115008/504984375 
8 407 30842 1 /3 1 6572300000 
8 -3209/4406250 
s -68947/9075000 
8 1053709375/1315607832 
8 664346654667/1317136718750 
8 126594108/1053709375 
8 "!J 9 S 4953001 , 6 4/62563994 140625 
- "12781897851723/50051195312500 

* If 7 53 2f 3259 / 73171 367 1 87500 

f 1°75 3 1 0068 14/658568359375 

?L TBc * G*( 5 ) » 0 * 101873 ID -03 

B* (-0.4512000D+01, 0) 
l! (0, 0.3606560D+01) 


STABILITY LIMITS 


TABLE Vr- Continued 
(b) Continued 


SIQMA a 0.15D0} A(6) s 0.60D0 


C(0) = 315161/1600000 

C(1) 80 

C(2) s 317696/16921075 

C(3) * 1 30^92289/362 1 55200000 

C(4) s -351/2000000 

C(5) 8 -305991/45100000 

C(6) s 1334375/1687314 

B(6,0) 8 103235259163/213500000000 

B(6,1) s 915471/13343750 

B(6,2) s -47130754/595703125 

B(6,3) s -914201623309/2348500000000 

B(6,4) s 75110391747/266875000000 

B(6,5) a 34464232107/146781250000 

SCALED NORM OP TEC, Q*(5) a 0. 30533 16D-03 

SCALED STABILITY LIMITS R; (-0.3734667D+01 , 0) 

I: (0, O. 36 O 656 OD+OI) 


SIQMA s 0.20D0J A(6) s 0.60D0 


C(0) s 801613/4050000 
C(1) s 0 

C(2) s 237568/24046875 

C(3) s -222404507/73641150000 

C(4) s 7487/4125000 

C(5) s 10624/1959375 

C(6) s 9284375/11782584 

B(6,0) s 22313787003/46421875000 

B(6,1) s -2651112/46421875 

B(6,2) s 15744308256/551259765625 

B(6,3) 8 -461266147641/882015625000 

B(6,4) s 10844473266/29013671875 

B(6,5) s 1723827744/5802734375 

SCALED NORM OP TEC, Q*(5) a 0.6400000D-03 

SCALED STABILITY LIMITS R: (-0.3273000D+01, 0) 

I: (0, 0.3606560D+01) 
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TABLE V,» Continued 
(b) Continued 


SIQMA b 0.25D0} A(6) ■ Q.6QD0 


C(0) s 40571/207360 

C(1) 8 0 

G(2) a -2656/115425 

C(3) ■ -27104389/2903662080 

0(4) b 889/163200 

' ') a 1279/36960 

Lv5) a 3086875/3874878 

0(6,0) s 1952712303/3951200000 

B<6,1) b -1542231/6173750 

0(6,2) s 292545258/1466265625 

0(6,3) = -536606220891/825800800000 

0(6,4) b 2267080407/4939000000 

0(6,5) * 942241167/2716430000 

SCALED NORM OP TEC, Q*(5) a 0, 1 101072D-02 

SCALED STABILITY LIMITS R} (-0.2962400D+01, 0) 

Is (0, 0.3606560D+01) 


SIGMA s 0.30D0; A(6) a 0.60DO 


C(0) s 9543/50000 
C(1) 8 0 

C(2) 8 -1078528/11578125 

C(3) a -93194543/5047350000 

C(4) r 54459/5125000 

C(5) 8 189243/2200000 

C(6) a 13578125/16476096 

B(6,0) s 4412352476/8486328125 

0(6,1) s -33320484/67890625 

0(6,2) s 338943067216/806201171875 

0(6,3) = -5382504160943/7094570312500 

0(6,4) s 89887658679/169726562500 

0(6,5) 8 35427725424/93349609375 

SCALED NORM OP TEC, G»(5) s 0. 1670540D-02 

SCALED STABILITY LIMITS R: (-0.2737667D+01 , 0) 

Is (0, 0.3606560D+01) 


TABLE V.« Continued 


(b) Continued 


3IQHA b Q.35DQ; A(6) s Q.4GD0 


C(0) a 830677/7200000 

C(1) a 0 

C(2) a 1 4337008/111859375 

C( 3) b 6313094879/765943200000 

C(4) b -143423/32250000 

C(5) a -5274409/59400000 

C(6) b 1114759375/1777574592 

B(6,0) a 35745803245261/150492515625000 

B(6, 1) a 369196744/1114759375 

B(6,2) a -962010407052128/1787098623046875 

B(6,3) * -21076530101281837/31452935765625000 

B(6,4) a 1712370111046/3483623046875 

B(6,5) * 4200945736554/7663970703125 

SCALED NORM OP TEC, G*(5) = 0.2290976D-02 

SCALED STABILITY LIMITS Rs (-0.25674290+01,0) 

Is (0, 0 . 3606560D+0 1 ) 


SIGMA e O.40D0; A(6) 8 0.40D0 


C(0) a 27781/225000 
C(1) a 0 

C(2) s 4030464/12765625 

C(3) 8 -4826809/2915550000 

C(4) 8 1114/984375 

C(5) s -40288/1753125 

C(6) s 13346875/22842288 

B(6,0) a 76210981292/225228515625 

B(6,1) 8 -11195344/66734375 

B(6,2) a -3146683828928/21396708984375 

B(6,3) * -42194535738914/47072759765625 

B(6,4) 8 26701738896/41708984375 

B(6,5) * 58049867904/91759765625 

SCALED NORM OF TEC, G*(5) 8 0.2880000D-02 

SCALED STABILITY LIMITS Rs (-0.2434750D+01 , 0) 

It (0, 0.3606560D+01 ) 
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TABLE V." Continued 
(b) Continued 


9IGMA r 0.45DQJ A(6) * 0.40PQ 


C(0) a 631 3/00000 
C(1) ■ 0 

C(2> a 939984/3899375 

C(3) a -79085409/6729800000 

C<4) a 376777/41000000 

C(S) a 1288629/17600000 

C(6) a 24690625/43936356 

B<6,0) a 7205960559/15431640625 

B(6, 1 ) a -97364088/123453125 

B(6,2) a 531961896544/1466005059375 

B(6,3) = -7208130979231/6450425781250 

B(6,4) a 120503849193/154316406250 

B(6,5) a 118047108816/169748046875 

SCALED NORM OP TEC, Q»(5) a 0.3502154D-02 

SCALED STABILITY LIMITS Rs (-0.2328889D+01 , 0) 

Is (0, 0.3606560D+01) 


SIQMA a O.50DOJ A(6) a 0.40D0 


C(0) a 721/5760 

C(1) a 0 

C(2) a 512/4275 

C(3) a -371293/18860160 

C(4) a 107/9600 

C(5) a 98/495 

C(6) a 30625/54144 

B(6,0) a 9933211/16537500 

B(6,1) a -43056/30625 

B(6,2) a 51378752/56109375 

B(6,3) a -4458159667/3456337500 

B(6,4) a 677552/765625 

B(6,5) » 1213728/1684375 

SCALED NORM OF TEC, G«(5) a 0.4207160D-02 

SCALED STABILITY LIMITS Rs (-0.2243600D+01 , 0) 

I: (0, 0.3606560D+01) 


TABLE V.- Continued 
(b) Continued 


SIGMA e 0.5500} A(6) s Q.40DQ 


C(0) 8 874507/7200000 

C(1) so 

C(2) s -5155568/82828125 

C(3) 8 -159284697/6946400000 

C(4) s 763147/58500000 

C<5) b 1459733/4200000 

C(6) s 447859375/742577472 

B(6,0) 8 42737353009567/60461015625000 

8(6,1 ) a -4377093688/2239296875 

8(6,2) s 1001659717578784/717974560546875 

8(6,3) a -1571898291387749/1148759296875000 

8(6,4) a 1293377129187/1399560546875 

8(6,5) s 195073554858/279912109375 

SCALED NORM OP TEC, 0*(5) s 0.5955773D-02 

SCALED STABILITY LIMITS R: (-0.2174364D+01 , 0) 

Is ( 0 , 0.3606560D+01) 


SIGMA a 0.60DO; A(6) s 0.10D0 


C(0) s 23929/150000 

C(1) = 0 

C(2) s 565248/296875 

C(3) a -6311981/19950000 

C(4) a 19899/125000 

C(5) a -8208/3125 

C(6) s 289/168 

8(6,0) s -1213233/1445000 

B(6,1) s 24/5 

8(6,2) s -61850336/17159375 

8(6,3) a 265378841/302005000 

8(6,4) a -437346/903125 

8(6,5) a 495936/1986875 

SCALED NORM OP TEC, G»(5) = 0.6109403D-02 

SCALED STABILITY LIMITS R{ (-0.21 17500D+01 , 0) 

It (0, 0.3606560D+01) 
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TABLE V.- Continued 
(b) Continued 


SIGMA u O. 65 DOJ A ( 6 ) b 0 . 1 QDQ 


C(0) a 2150839/14400000 

C(1) B 0 

0(2) b 47449792/29390625 

C(3) b -27797593031/71227200000 

0(4) a 7951619/42000000 

0(5) b -16920449/9900000 

0(6) b 112799/98406 

B(6,0) a -366959163251/487291680000 

B(6, 1 ) 8 4844333/1127990 

B(6,2) a -588204877546/180830896875 

B(6,3) * 6048781368877/9258541920000 

B(6,4) s -725826801/22559800000 

B(6,5) 8 428189853/1127990000 

SCALED NORM OP TEC, G«(5) s 0.5947044D-02 

SCALED STABILITY LIMITS R: (-0.2070923D+01 , 0) 

I: (0, 0.3606560D+01) 


SIGMA s O.70DO; A(6) 8 0.10D0 


C(0) s 126439/900000 

C(1) s 0 

C(2) s 5406464/3859375 

C(3) * -2747196907/5454900000 

C(4) s 65219/281250 

C(5) 8 -1016407/825000 

C(6) s 26123/27144 

B(6,0) s -6659987/11677550 

B(6,1) s 441028/130615 

B(6,2) s - 1 05258859888/4 078434375 

B(6,3) * 74987588767/184265111250 

B(6,4) 8 -50534211/326537500 

B(6,5) 8 81719274/179595625 

SCALED NORM OP TEC, 0*(5) a 0.5457097D-02 

SCALED STABILITY LIMITS R: ( -0. 20327 14D+01, 0) 

I: (0, 0.3606560D+01) 


* 1 1' ‘‘ 


r;r t\yt 


TABLE V.- Continued 
(b) Continued 


SIGMA s 0.75D0; A(6) s O.IODO 


C(0) c 1019/7680 
C(1> 8 0 

C(2) s 352/285 

C(3) * -371293/53501*0 

C(4) s 189/640 

C(5) = -1647/1760 

C(6) a 29/30 

B(6,0) s -2463/7424 
B(6, 1 ) s 129/58 

B<6,2) s -4206/2755 
B(6,3) * 327691/1551616 

B(6,4) s -279/9280 
B(6,5) * 2313/5104 

SCALED NORM OF TEC, G«(5) s 0.4661 104D-02 

SCALED STABILITY LIMITS R: (-0.20012000+01 , 0) 

I: (0, 0.3606560D+01) 


SIGMA 8 0.80D0; A(6) = O.IODO 


C(0) 8 393/3125 

C(1) s 0 

C(2) s 50069504/45421875 

C(3) s -24876631/23512500 

C(4) s 18544/46875 

C(5) 8 -227072/309375 

C(6) s 715/612 

B(6,0) s -650311/12065625 

B(6,1) s 3272/3575 

B(6,2) s -415033376/1146234375 

B(6,3) = 18109624/193978125 

B(6*4) s 76432/2234375 

B(6,5) * 1834368/4915625 

SCALED NORM OF TEC, G*(5) 8 0. 36203870-02 

SCALED STABILITY LIMITS R: (-0.19747500+01,0) 

I: (0.1676275D+0, 0.2094762D+01) 


TABLE V.- Continued 
(b) Continued 


SIOMA a 0.85DO; A(6) ■ 0.10D0 ' 


C(0) s 1721 699/ 14400000 

C(1) * r 

0(2) s 501157088/50765625 

C(3) b -3723113759 19 /1 9060800000 
C(4) s 10236091/18000000 

C(5) a -1944103/3300000 
C(6) s 12089/6498 

8(6,0) 8 14522782003/52224480000 
B(6, 1 ) s -729 13/120890 
B(6,2) s 146868746/145715625 
B(6,3) 8 24076851851/574469280000 
B(6,4) 8 104040561/2417800000 
B(6,5) s 308770881/1329790000 
SCALED NORM OF TEC, G # (5) * 0.2442667D-02 

SCALED STABILITY LIMITS R: (-0.1952000D+01, 0) 

(0.1501106D+0, 0.2144200D+01) 


T* 


SIOMA s 0.90D0; A(6) s 0.10D0 


C(0) s 2849/25000 
C(1) a 0 

C(2) * 1879296/2078125 

C(3) b -4455516/653125 
C(4) s 114939/125000 
0(5) a -265113/550000 
0 ( 6 ) s 713/112 
B(6,0) s 629873/891250 
B(6,1) s -8892/3565 
B(6 ,2) s 113733424/42334375 
B(6,3) s 214123/12017500 
B(6,4) 8 144081/8912500 
B(6,5) s 328356/4901875 

!«*}•“ “ORM OF TEC. 0»(5) • 0.1268702D-02 

SCALED STABILm LIMITS Rt 0) 

l! (0.1371556D+0, 0.2158556D+01) 


TABLE V.~ Continued 


(b) Concluded 


SIGMA = 0.95DO} A(6) s 0.1D0 


C(0) s 1567469/14400000 

C(1) 8 0 

C(2) a 894368/1078125 

C(3) a 7146276371/1108800000 

C(4) s 11800007/6000000 

C(5) = -11919137/29700000 

C(6) s -69107/8694 

B(6,0) s 384031234573/298542240000 

B(6,1) s -3424921/691070 

B(6,2) s 27848646422/5830903125 

B(6,3) a -128214678219/3283964640000 

B(6,4) a -201293489/13821400000 

B(6,5) a -401374449/7601770000 

SCALED NORM OF TEC, G»(5) s 0.3789098D-03 

SCALED STABILITY LIMITS R: (-0.19104210*01, 0) 

It (0.12731680*0, 0.21576530+01) 





TABLE V,- Continued 
(o) Scaled, fifth-order algorithm 


C^k.X 


for X a where l 

are the unaoaled ooeffiolenta (table V(a)) K,A 


SIGMA a 0.5000000D-01 { 


A( 1 ) : 

A(2) « 

A(3) : 

A(4) : 

A(5) = 

A(6) : 

A(7) * 

C(0) a 
C(1) s 

0 ( 2 ) a 

C(3> a 
C(4) 8 

C(5) a 
C(6) a 
C(7) a 

B0(6) a 
B(6, 1) 8 
B(6,2) a 
B(6,3) a 
B(6,4) a 
B(6,5) a 
B0(7) a 
B(7,1) a 

fi(7,2) a 

B(7,3) a 
B(7,4) a 
B(7,5) a 
B(7,6) a 
SCALED 
SCALED 


> . 5000000000000000D+0 1 
» • 7500000000000000D+0 1 

! *18461538461538460+02 
.2000000000000000D+02 
. 1 0000000000000000+02 
. 3500000000000000D+00 
. 8500000000000000D+00 
. 1049397167478 102D+00 
. OOOOOOOOOOOOOOOOD+OO 
-.3224643108141 323D-04 
-. 1766214053432440D-05 
. 94240 1 3548044395D-06 
. 1535439210883020D-04 
.52142847184935 37D+00 
. 3736495272545072D+00 
.336181 0372988544D+00 
.13673071478920280-01 
. OOOOOOOOOOOOOOOOD+OO 
.65215681 95420585D-02 
-.63756769731952250-02 
. OOOOOOOOOOOOOOOOD+OO 
- . 3 65858087747200 1 D +00 
- . 201453450608 1 222D-0 1 
. 1 08764498610 1 98 1 D -02 
-.12816921096675150-01 
.11381 92692070809D-0 1 
.311 690752 1 207 1 85D-02 
.12332338744766700+01 

0F TEC * Q # ( 6 ) a .25466970-06 

STABILITY LIMITS R: (-.1114D+02, 0.0) 

Is (0.0, 0.0) 
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TABLE V.- Continued 


(o) Continued 


SIQMA s 0. 10000000+00 


A( 1 ) b . 25000000000000000+0 1 
A(2) s . 3750000000000000 D+0 1 
A(3) s .9230769230769231D+01 
A ( 4 ) 8 .1 000000000 OOOOOOD+02 

A(5) 8 .50000000000000000+01 

A(6) 8 .25000000000000000+00 

A(7) 8 .80000000000000000+00 

C(0) 8 .44512744708994700-01 

C(1) s .00000000000000000+00 
C(2) s -.38095363239953500-03 
C( 3) « -.26115732478398960-04 

C(4) s .14035886818495500-04 
C(5) s .20463453004806370-03 
C(6) s .48006256452425410+00 
C(7) * .47561308971476260+00 

BO (6) s .22094115484340500+00 
B(6, 1 ) s .41334586699043540-01 
B(6,2) s .00000000000000000+00 
B(6,3) * -.6587223327993963D-01 
B(6,4) s .53596491737491080-01 
B(6,5) 8 .00000000000000000+00 

BO (7) s -.42956730219763850+00 
B(7» 1 ) 8 -.47512231125677210-01 
B(7,2) 8 .56512311382852310-02 

B(7 ,3) * .46750866966850420-01 

B(7,4) 8 -.40845781331627700-01 
B(7,5) s .1644610650843368D-01 
B(7*6) s .12490771100413740+01 
SCALED MORM OP TEC, Q»(6) a .4513416D-05 

SCALED STABILITY LIMITS R: (-.70300+01, 0.0) 

I* (0.0, .36040+01) 


SIQMA 8 0.15000000+00 


A( 1 ) & .16666666666666670+01 
A(2) s . 2500000000000000D+0 1 
A(3) s .61538461538461540+01 
A(4) 8 .66666666666666670+01 
A(5) 8 .33333333333333330+01 
A(6) a .45000000000000000+00 
A(7) 8 .90000000000000000+00 



TABLE V.~ Continued 
(q) Continued 


SIGMA o 0.15QQ000D+0Q = Conoluded 


C(0) b .14015950500540860+00 
C( 1 ) B .00000000000000000+00 
C(2) 8 .49865952686594560-02 

C(3> b .23808821004188780-03 
C(4) a -.12627479007223780-03 
C(5) s -.22204764400220900-02 
C(6) b .59941125645508540+00 
C(7) 8 .24955930549089900+00 

BO (6) s .51202713255560840+00 
B(6,1 ) a -. 169924328 10 16840D+00 

B(6,2) s .00000000000000000+00 

B(6 , 3 ) 8 .65297197576850930+00 
B(6,4) 8 -.54507470022243370+00 
B(6,5) a .00000000000000000+00 
B0(7) s -.72521205847745730+00 
B(7,1) 8 .48359165576957000+00 
B(7,2) s -. 791760429587729 3D-01 
B(7»3) a -.15247093839287420+01 
B(7»4) 8 .1279935176012341D+01 
B(7,5) a -.43792977077363620-02 
B(7»6) a .1 46994995 1290798D+01 
SCALED NORM OF TEC G»(6) 8 .11625570-03 

SCALED STABILITY LIMITS Rj (-.43330+01,0.0) 

I: (0.0, 0.0) 


SIGMA s 0.20000000+00 


A(1 ) 
A(2) 
A(3) 
A(4) 
A(5) 
A(6) 
A(7) 
C(0) 
C(1) 
0 ( 2 ) 
C(3) 
C(4) 
C(5) 
C(6) 


a .12500000000000000+01 
a . 1 8750000000000000+0 1 
a .46153846153846150+01 
a .50000000000000000+01 
a .25000000000000000+01 
a .41666666666666670+00 
8 . 1 0000000000000000+0 1 
8 .12498952380952380+00 
8 .00000000000000000+00 
8 -. 17412801 1 1844235D-01 
8 -.17266412915609770-03 
8 .82891774891774800-04 
a .4130724386724384D-02 
a .62687602569482240+00 
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TABLE V.- Continued 
(o) Continued 


SIGMA a 0.20QQQQQD+QQ - Conoluded 


C(7) s .26 1506299581 636 1D+00 
B0(6) s .32545914297783460+00 
B(6,1) s .10255254655011360+00 
B(6,2) b , OOOOOOOOQOOOOOOOD+OO 
B(6,3) s -.39883965353235320-01 
B(6,4) s .2853894249195378D-01 
B(6,5) s .0000000000000000D+00 
BO (7) = -.242521 1742757191D+00 

B(7,1) * -.30829783203501400+00 
B(7 ,2) 8 .10871822499691550+00 

B(7,3) * .57266104666919630-02 
B(7»4) s -.90433530513206510-02 
B(7,5) 8 .47066194541054890-01 
B(7,6) 8 .13983513293573910+01 

SCALED NORM OF TEC, G»(6) a .3767397D-04 

SCALED STABILITY LIMITS Rt (-.59100+01, 0.0) 

I: (0.0, 0.0) 


SIGMA a 0.25000000+00 


A ( 1 ) a .1 OOOOOOOOOOOOOOOD+O 1 
A(2) a .15000000000000000+01 
A(3) s .3692307692307692D+01 
A(4) s .40000000000000000+01 
A(5) s .2000000000000000D+01 
A(6) s . 2500000000000000D+00 
A(7) s .10000000000000000+01 
C(0) s 6557126322751 310-01 
C( 1 ) 8 .00000000000000000+00 
C(2) 8 -.20930548593706490+00 
C(3) * -.3710393274856326D-02 
C(4) s .18991402116402110-02 
C(5) s .56018089053803330-01 
C(6) s .57402149338941040+00 
C(7) s .56452003023431590+00 
B0(6) s -.12038596976257230+00 
B(6, 1 ) b .6481601 383594973D+00 
B(6,2) a .OOOOOOOOOOOOOOOOD+OO 
B(6,3) 8 -.16061062420916590+01 
B(6,4) s .13283320734947340+01 
B(6,5) s .OOOOOOOOOOOOOOOOD+OO 


TABLE V.- Continued 
(a) Continued 


SIGMA S G.250QGQQB+QQ - Conoluded 


00(7) 
0(7,1) 
B(7,2) 
0(7,3) 
0(7,10 
0(7,5) b 

0(7,6) s 
SCALED 
SCALED 


8 - 


8 - 


*♦2 1 3556639011 306 1 0+00 
10143132581742230+01 
53919947671007820+00 
12930168081 17887D+01 
11277968667H9287D+01 
.12591489094950250+00 
.762623285211173560+00 
NORM OP TEC, G*(6) s .161126040-02 

STABILITY LIMITS R: (-.25000+01,0.0) 

li (0.0, .25160+01) 


SIGMA & 0.30000000+00 


A( 1 ) 
A(2) 
A(3) 
A(4) 
A(5) 
A(6) 
A(7) 
C(0) 
C(1) 
C(2) 
C(3) 
C(4) 
C(5) 
C(6) 
C(7) 
00 ( 6 ) 
0 ( 6 , 1 ) 
0 ( 6 , 2 ) 
0(6,3) 
0(6,4) 
0(6,5) 
00(7) 
0(7,1) 
0(7,2) 
0(7,3) 
B(7 , 4 ) 
0(7,5) 


8 .83333333333333330+00 
8 . 1 2500000000000000+0 1 
b .30769230769230770+01 
8 * 3333333333333333O+0 1 
b .16666666666666670+01 

* .37500000000000000+00 
b . 1 OOOOOOOOOOOOOOOD+O 1 

* .11099412698412700+00 
8 . 00000000000000000+00 
8 -.16137037450769780+00 
s -.27685847105151570-03 

* .11287381431445740-03 
8 .17750096355257630-01 

* .57629605671192350+00 
8 .45649407911312670+00 
a .27154042614440980+00 
= .111 96000647879360+00 

b .ooooOoooooooooooo+oo 

8 -.20854353045564240-01 
8 .12353920422360910-01 
8 .00000000000000000+00 
8 .57514950353042260-01 

* -.69192112924813260-01 

* .17988263604295070+00 
8 -.11713976386916620+00 
8 .78501841048195260-01 
8 .81407899141567190-01 


TABLE V.- Continued 
(g) Continued 


SIQMA r 0 . 3Q0QQ0QD+0Q - Conoluded 


B(7,G) » . 7 89024550208224 1 0+00 

SCALED NORM OF TEC, 0»(6) a .23477860-03 

SCALED STABILITY LIMITS R: (-.40400+01, 0.0) 

I: (.32130+01, .40140+01) 


SIOMA a 0. 3500000D4-00 


A( 1 ) a .71428571428571420+00 
A(2) a .10714285714285710401 
A(3) * .2637362637362637 D401 

A(4) S .28571428571428570401 
A(5) a . 1428571 42857 1429D401 
A(6) a .35000000000000000400 
A(7) a ,10000000000000000+01 
C(0) a .10322866732804240400 

C(1) a .00000000000000000400 
C(2) a -.8356817496767217D400 
C(3) * -.31897229734589710-03 

C(4) a .11635039082474970-03 
C(5) a .34439274198943050-01 
C(6) a .54082042160191180400 
C(7) a .1157396008454346D401 
B0(6) a .25072341721558560400 
B(6, 1 ) a .10333576725023070400 
B(6,2) a .00000000000000000400 
B(6,3) = .43843455076373380-02 

B(6,4) a -.84435299734536790-0? 

B(6,5) a ,00000000000000000400 
B0(7) a .19552008974590990400 
B(7,1) a .36160553108183600400 
B(7 ,2) a .11802002527466480+00 
B(7,3) a -.90345086345623350-01 
B(7,4) a .61588147668645840-01 
B(7,5) a .49883552480922830-01 
B(7,6) a .30372774009364400+00 
SCALED NORM OF TEC 0*(6) a .49772490-03 

SCALED STABILITY LIMITS R: (-.35140+01,0.0) 

I? (.28140+01, .34920+01) 


TABLE! V.- Continued 
(o) Continued 


9IQMA r 0 . 40000000+00 


A(1) ■ .62500000000000000+00 

A(2) a . 9 3750000000000000+00 
A(3) b .23076923076923080+01 
A(4) a . 2500000000000000D<<<01 
A(5) a .12500000000000000+01 
A(6) a . 3 1 25000000000000D-<-00 
A(7) s . 1000000000000000D+01 
C(0) a .90988359788359730-01 
CO) 8 .00000000000000000+00 
C(2) a .151 16897666388190401 
C(3) a -.12332042111951820-02 
C(4) a .55176417233559870-03 
C(5) s .92189552669552560-0 1 
C(6) a .48988132736263550+00 
C(7) a -.11840675664205070+01 
BO (6) a .19974879950871270+00 
B(6,1) a .13499306819994130+00 

B(6,2) a .00000000000000000+00 

B(6,3) a -.10432305856269220+00 
B(6,4) a .82081190854038210-01 
B(6,5) s .00000000000000000+00 
B0<7) a .24002487612026710+00 
B(7 , 1 ) s .13420294634216280+01 
B(7,2) a -.26922744501219690+00 
B<7,3) * .12122297761278200+00 
6(7*4) a -.72168790140206940-01 
B(7 ,5) a -.77443418043657470-01 
B(7,6) 8 -.28443766395861560+00 
SCALED NORM OP TEC, G*(6) a .85754430-03 

SCALED STABILITY LIMITS R: (-.34250+01, 0.0) 

I: (.25860+01, .35810+01) 


SIOMA a 0.45000000+00 


A(D a .55555555555555560+00 

A<2) a .83333333333333330+00 

A(3) s .20512820512820510+01 

A(4) a .2222222222222222 

A(5) a .11111111111111110+01 

A(6) s .27272727272727270+00 

A(7) 8 .10000000000000000+01 


TABLE V.- Continued 
(c) Continued 


SIGMA b Q.4GQQQQQD+QQ - Concluded 


CJ(O) • .77561035317460980-01 
C( 1 ) ■ . QQQQOQQOQQQGQQQQB+O0 

C(?,) ■ .802861 363 137S733P+00 

0(3) ■ - . 292289007 178890 1 0-02 

C(M a .13581003633671990-02 
0(5) b .27682950868*108670+00 
0 ( 6 ) b . A 372 1 409925 1 80530+00 
0(7) b -.59290201668250450+00 
BO ( 6 ) b .16032840689852820+00 
B(6,1) b .14448738836011110+00 

B(6,2) s .00000000000000000+00 

8(6,3) s -.16512876420847860+00 
B(6,4) a .13304024167711200+00 
B ( 6 , 5 ) 8 . 00000000000000000+00 
BO (7) a -.58134168340624750-01 
B(7 , 1 ) a .30231910157380380+01 
B(7 , 2 ) a -.14334505950137020+01 
B(7,3) s .56208833311652480+00 
B(7,4) b -. 335177927 1591276D+00 
B(7 ,5) 8 .—.22221574956004930+00 
B(7,6) s -.53630090878105890+00 
SCALED NORM OP TEC, G«( 6 ) 8 .15008750-02 

SCALED STABILITY LIMITS R: (-.32020+01,0.0) 

I: (.23600+01, .34070+01) 


SIQMA s 0.5000000D+00 


A(1) s .50000000000000000+00 
A(2) s .75000000000000000+00 
6(3) a .18461538461538460+01 
A(4) s .20000000000000000+01 
A ( 5 ) 8 .1 0000000000000000+0 1 
A(6) 8 .52000000000000000+00 
A ( 7 ) s . 47000000000000000+00 
C(0) s .10988545536417930+00 
C(1) 8 .00000000000000000+00 
C<2) 8 .10441584983853700+01 
C(3) s .79654166757858620-02 
C(4) a -. 40482078717372860-02 
C(5) s -.65788273571292570-01 
C(6) s -.36123662256959460+01 
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TABLE V,« Continued 
(a) Continued 


SIQMA » O.PQQQQQQD+OO - Conoluded 


cm - .35201933367136410+01 

00(6) b .233697366445030 1B+0Q 
0(6,1) b .20921240548464330+00 
0(6,2) ■ .00000000000000000+00 
0(6,3) b .23312502739335920-01 
B(6,4) ■ -.26222354669009380-01 
0(6,5/ b .00000000000000000+00 
00(7) s .21856275582570970+00 
0(7,1) b .20112318445008030+00 
0(7,2) b -.50173732613543270-01 
0(7,3) b .50579049669133520-01 
0(7,4) a -.43365855281240700-01 
0(7,5) b -.15893917120429160-01 
0(7,6) s .10916851507028960+00 
SCALED NORM OF TEC, Q«(6) a .32957250-03 

SCALED STABILITY LIMITS Rs (-.34020+01,0.0) 

I: (0.0, .11810+01) 


SIOMA b 0.55000000+00 


A( 1 ) s .45454545454545450+00 
A(2) s .68181818181818180+00 
A(3) = .16783216783216780+01 
A(4) b .18181818 181818 180+01 
A(5) s .90909090909090910+00 
A(6) s .26000000000000000+00 
A(7) s .83000000000000000+00 

c(o) s .77850291891007090-01 

C(1) 8 .00000000000000000+00 
c( 2 ) 8 .39793740138230730+00 
C(3) b .20191174553617090-02 
C(4) s -.10286559790783990-02 
C(5) b .22130413561593210+00 
C(6) b .39406434616719190+00 
C(7) b -.92146636532721630-01 
B0(6) b .13049315899263240+00 
0(6,1) b .16950419790553560+00 
B(6,2) s .000000000000000004-00 
0(6,3) b -.21074699648319390+00 
0(6,4) b .17074963958502580+00 
0(6,5) a .OOOOOOOOOOOOOOOOD+OO 



' It M '{ * 


1 

4 

TABLE V.- Continued | 

. U) . . c °" “! uad 4 

SIGMA in 0.5500000D+00 - Conoluded I 


B0(7) b -.12186124450129080+01 
B(7»1 ) a , 1 1 698309445846960+02 
B<7,2) b -.604657764 1 630584D+0 1 
B(7,3) b .30632539455940240+01 
B(7,4) b -.18461568989224580+01 
B(7,5) a -.21283877523525450+01 
B(7,6) b -.26918286535224870+01 
SCALED NORM OP TEC, G»(6) s .3221258D-02 

SCALED STABILITY LIMITS R{ (-.2764D+01, 0.0) 

I: (.1988D+01, .29910+01) 


SIGMA b 0.60000000+00 


A(1) s .4166666666666667D+00 
A(2) s .62500000000000000+00 
A(3) 8 .15384615384615380+01 
A(4) 8 .16666666666666670+01 
A(5) 8 .8333333333333333D+00 
A(6) 8 .33000000000000000+00 
A(7) s .7200000000000000D+00 
C(0) 8 .93560344717011490-01 
C(1) s .OOOOOOOOOOOOOOOOD+OO 
C(2) s -.64892741212531840+00 
C(3) s .14821798325568710-01 
C(4) s -.71136124276832990-02 
C(5) 8 -.24774011808216770-02 
C(6) s .56256598953145810+00 
C(7) 8 .98757029315978520+00 
BO(6) 8 .13826261139150720+00 
B(6,1 ) s .23469110602695830+00 
B(6,2) s .00000000000000000+00 
B(6,3) * -.22036223185243670+00 
B(6,4) s .17740851443397130+00 
B(6,5) a .00000000000000000+00 
B0(7) 8 .16914919246580490+00 
B(7,1) 8 .16861776577261570+00 
B(7<2) 8 -.27257766741746070-02 
B(7,3) s -.529653224936337 3B-01 
8(7,4) a .2281937886921284D-01 
B(7,5) a .19294261579362270+00 




A(1) a .38461538461538460+00 
A(2) a .5769230769230769D+00 
A(3) * .14201183431952660+01 

A(4) 8 .15384615384615380+01 

A(5) * .76923076923076920+00 

A(6) 8 .450000000000000 0D+00 

A(7) = .60000000000000000+00 

C(0) a ,1019736151 528440D+00 
C(1) a .00000000000000000+00 
C(2) a -.16646464637823530+02 
C(3) * .35345607520267780-01 

C( 4 ) a - . 1 6550629702 1 1 5460-0 1 
C(5) a -.16978109903605060+00 
C(6) a .27179218702222960+01 
C(7) = .14977555273666290+02 

BO ( 6 ) a .1 62034688860 1 6020+00 
B(6,1) a .301 273283227 1356D+00 
B(6,2) a .00000000000000000+00 
B(6,3) a -.49427966646655030-01 
B(6,4) a .36119994559359200-01 
B(6,5) a .00000000000000000+00 
80(7) a .13595038554913840+00 
B(7,1) a . 4594 121 042544327D+00 
B(7,2) a -.23957460722224670-01 
B(7, 3) * .13253632436263990-01 
B(7»4) a -.87615451351783380-02 
B(7,5) a .2931812700755386D-02 
B(7»6) a .21171070916812470-01 
SCALED NORM OP TEC, G*<6) a .17751570-02 

SCALED STABILITY LIMITS R; (-.32090+01,0.0) 

I: (.21080+01, .36360+01; 










TABLE V.» Continued 
(o) Continued 


SIGMA b 0.7000000D+00 


A(1) b .35714285714285710+00 
A(2) s .53571428571428570+00 
A(3) a .13186813186813190+01 
A(4) a . 1 42857 1 42857 1 429 D+0 1 
A(5) s .71428571428571430+00 
A(6) s . 4600000000000000D+00 
A(7) s .59000000000000000+00 
C(0) a .993861927181381 OD-O 1 
C(1) s .0000000000000000D+00 
C(2) s -.1164829345465175D+02 
C(3) a . 56069477969 13903D-01 
C(4) s -.2523500352943741D-01 
C(5) s -.30703678176574100+00 
C(6) s .4963332454942884D+01 
C(7) s .78617771143167680+01 
B0(6) s .1571494738728330D+00 
B<6,1) s .3031440894261 848D+00 
B(6,2) s .OOOOOOOOOOOOOOOOD+OO 
B(6,3) * .1862 19677478696 ID-01 

B(6,4) s -.18915531046887390-01 
B(6,5) s .OOOOOOOOOOOOOOOOD+OO 
B0(7) s .1009320434329384D+00 
B(7,1) 8 . 5 1 2207 58704744 1 8D+00 

B(7,2) * -.78027276554405780-01 
B(7»3) 8 .15374880434390030-01 

B(7 *4 ) 8 -.44534483398697800-02 
B(7,5) s -.47895389861884590-02 
B(7 f 6) a .4875575296569375D-01 
SCALED NORM OF TEC, G«(6) s .14301450-02 

SCALED STABILITY LIMITS R: (-.39070+01,0.0) 

I: (0.0, 0.0) 


SIGMA s 0.75000000+00 


A(1) 

S 

.33333333333333330+00 

A(2) 

s 

. 50000000000000000+00 

A(3) 

s 

.12307692307692310+01 

A(4) 

s 

.13333333333333330+01 

A(5) 

s 

.66666666666666670+00 

A(6) 

wm 

. 45000000000000000+00 

A(7) 

s 

. 59000000000000000+00 
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TABLE V.- Continued 


(q) Continued 


SIGMA b 0.7500000D+00 - Conoluded 


C<0) a .958763858061 99 13D-01 
C(1) a .00000000000000000+00 
C(2) a -„ 993425786688945 10+01 
C(3) a .8389383603325926D-01 
C(4) a 35929 10703709527D-01 
C(5) a -.5843105264300945D+00 
C(6) a .64983448994177600+01 
C(7) a .48763823790994220+01 
80(6 ) a .15052704756036770+00 
B(6,1) a .29216839874675620+00 
B(6,2) a .OOOOOOOOOOOOOOOOD+OO 
B(6,3) a .57318993934346870-01 
B(6,4) a -.50014440241470180-01 
B(6,5) a .00000000000000000+00 
B0(7) a .34044537252968510-01 
B(7,1) * .6907242066214459D+00 

B(7,2) a -.2264502174892697D+00 
B(7,3) * .4487947377765493D-02 

B(7»4) a .1726173816542408D-01 
B(7,5) a -.13730199340895430-01 
B(7,6) a .83661987412561160-01 

SCALED NORM OF TEC, G*(6) a .36517270-02 

SCALED STABILITY LIMITS R: (-.52010+01,0.0) 

I: (0.0, 0.0) 


SIGMA a 0.80000000+00 


A(1) 

S 

. 3 1 250000000000000+00 

A(2) 

m 

. 46875000000000000+00 

A(3) 

s 

.11538461538461540+01 

A(4) 

8 

.12500000000000000+01 

A(5) 

8 

. 62500000000000000+00 

A(6) 

8 

. 4300000000000000D+00 

A(7) 

8 

. 58000000000000000+00 

C(0) 

8 

.91 7770 1 437947632D-0 1 

C(1) 

8 

.oooooooooOooooooo+oo 

C(2) 

8 

-.88999111581507100+01 

C( 3) 

8 

.11909993304736900+00 

C(4) 

8 

- . 47796500702064620-0 1 

C(5) 

8 

-. 767488385 1950577D+00 

C(6) 

8 

.67911256913504740+01 
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TABLE V.- Continued 
(a) Continued 


SIGMA B G.8Q00Q0QD+00 - Concluded 


C(7) « 

B0(6) = 

B(6 f 1 ) a 
B(6,2) * 
B(6*3) 8 
B(6,4) s 
B(6,5) * 
BO (7) s 
B(7,1) = 
B(7,2) ■ 
B(7»3) * 
B(7,4) 8 
B(7 *5) 8 
B(7,6) 8 
SCALED 
SCALED 


.37 131 9340527051 3D+01 
. 1 321587392059206D+00 
.296471 549 3885388D+00 
. OOOOOOOOOOOOOOOOD+OO 
.1985878378477856D-01 
-.18489072379237 

. OOOOOOOOOOOOOOOOD+OO 
- . 27998710660 1 3902D-02 
.81 9819 3282626641 D+00 
- . 37470 1 3763688526D+00 
.71 13607481 1031 13IMJ1 
- . 3002685456 127201 D-0 1 
. 807667 4027 364075D-02 
. 88496024895079 14D-01 
NORM OP TEC, 0*(6) * ^TOWID-OJ 

STABILITY LIMITS R: (-.2700D+01, 0.0) 

I: (.1606D+01, .3120D+01) 


SIGMA = 0.8500000D+00 


A( 1 ) 

8 

A(2) 

S 

A( 3) 

s 

A(4) 

s 

A(5) 

s 

A(6) 

s 

A(7) 

s 

C(0) 

s 

C(1) 

s 

C(2) 

s 

C(3) 

s 

C(4) 

8 

C(5) 

s 

C(6) 

s 

C(7) 

s 

B0(6) 

s 

B(6,1 ) 
B( 6,2) 

s 

s 

B(6,3) 

s 

B(6,4) 

s 

B(6,5) 

8 


.2941176470588235D+00 
.441 1764705882353D+00 
.1085972850678733D+01 
.1176470588235294D4-01 
. 588235294 1176471 D+OO 

. 4200000000000000D+00 
. 55OOOOOOOOOOOOOOD+OO 
.879 1034006448265D-01 
. OOOOOOOOOOOOOOOOD+OO 
-.1432598557850752D+O2 
. 1637724502289815D+00 
-.6025067094271 176D-01 
-.605122312396251 1D+00 
.1221191 469286895D+02 
. 352776 1 078684070D+0 1 
.1184034763470802D+00 
.30 12366421 08 1962D+00 
.OOOOOOOOOOOOOOOOD+OO 
.9087546933044339D-02 
-.8727665388320754D-02 
. OOOOOOOOOOOOOOOOD+OO 
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TABLE V.- Continued 
(o) Continued 


SIOMA 5 0.05000000+00 - Conoludod 


BO (7) a . 0053309076 199897D-02 
B(7,1) a .72275883007333000+00 
B(7,2) b -.31606326250504890+00 
B(7,3) a .53399835561240190-01 
B(7,4) a 19746867880022050-01 
B(7,5) s .60579463123815480-02 
B(7,6) s .94740209361903290-01 
SCALED NORM OF TEC, 0«(6) a .36119390-02 

SCALED STABILITY LIMITS R: (-*27880+01,0.0) 

I: ( . 1605D+01 , ,32320+01) 


SIGMA s 0.9000000D+00 


A( 1 ) = .27777777777777780+00 
A(2) s .41666666666666670+00 
A(3) a .10256410256410260+01 
A(4) s .11111111111111110+01 
A(5) s .55555555555555560+00 
A(6) s .41000000000000000+00 
A(7) 8 .50000000000000000+00 
C(0) s .83762758420442710-01 
C( 1 ) 8 .00000000000000000+00 
C(2) 8 -.46329178466165540+02 
C(3) s .21172843050175490+00 
C(4) a -.68158491695309560-01 
C(5) a .13725175770792540+00 
C(6) 8 .43290548093282350+02 
C(7) a .36740459179483780+01 
B0(6) s .10575517816709950+00 
B(6, 1 ) s .30502484906587520+00 
B(6,2) s .00000000000000000+00 
B(6,3) a -.21945945645767250-02 
B(6,4) 8 .141^5673316020210-02 
B(6 ,5) a .00000000000000000+00 
BO (7) s .97275839422458320-01 
B(7 , 1 ) s . 3^84586 1 060377080+00 
B(7 ,2) s -.74261310699345270-02 
B(7,3) a .3707095614794195D-01 
B(7,4) s -.24032246868840230-01 
B(7,5) s -.36106086630303260-01 
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TABLE V.- Conoluded 


(a) Cortaludod 


SIGMA a 0,90000000+00 - Conoluded 


B(7,6) = .121175905839496700+00 

SCALED NORM OF TEC, G«(6) s .3609105D-02 

SCALED STABILITY LIMITS R: (-.22690+01, 0.0) 

I: (0.0, . 1 142D+01 ) 


SIGMA s 0. 95000000+00 


A ( 1 ) s .2631 57894736842 1 D+00 

A(2) s .39473684210526320+00 

A(3) * .97165991902834010+00 

A(4) s .1052631578947368D+01 
A(5) s .52631578947368420+00 
A(6) a . 3750000000000000D+00 
A(7) s .45000000000000000+00 
C(0) s .78756132789580270-01 
C(1) s .00000000000000000+00 
C(2) s -.13918934655580860+02 

C(3) s .2491545555366749D+00 
C(4) s -.61966454571424580-01 
C(5) s .76777378408198430+00 
C(6) 8 .11133177189582230+02 

C(7) a .27520394481618090+01 
B0(6) s .11078585160391810+00 
B ( 6 , 1 ) s .1 600964963528354D+00 
B ( 6 , 2 ) = . OOOOOOOOOOOOOOOOD+OO 

B(6,3) s . 304837147089 1898D-01 
B(6 ,4) s -.26366062665672530-01 
B(6,5) s .OOOOOOOOOOOOOOOOD+OO 
B0(7) * .15709285674354190+00 

B(7, 1 ) * -.194939973928981 ID-01 
B(7,2) s .259296781 97843 16D+00 
B(7,3) s -.35633814808485140+00 
B(7,4) 8 .26036957964454320+00 

B(7,5) s .1 340898607825 139D+00 
B(7 ,6) «. .14983066328718960-01 

SCALED NORM OF TEC, G*(6) s .19106570-01 

SCALED STABILITY LIMITS R; (-.19460+01, 0.0) 

I: (.12600+01, .2208D+01) 



